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A COMPARISON OF SOLAR ATMOSPHERIC MODELS BY 
STARK EFFECT IN H-GAMMA : NITROGEN ABUNDANCE 


D. Mugglestone 
(Received 1959 September 7) 


Summary 


Three solar atmospheric models are considered, i.e., the models of Claas, 
Vitense and Swihart. The criterion of reliability adopted here is that the 
models should correctly predict the observed intensities in the far wing of 
the H-gamma absorption line. An improved theory of Stark broadening 
is used, due to Dr A. C. Kolb, where the ions are treated according to the 
usual Holtsmark statistical theory and also the effect of the electrons, and their 
finite velocities, is included. It is found that the electronic and ionic contri- 
butions to the absorption coefficient are of the same order of magnitude, in 
the far wing, and therefore the electronic effect should not be neglected. 
Results indicate that, of the models considered, the Swihart model is to be 
preferred. 

In a previous paper (g) the nitrogen abundance in the solar atmosphere 
was determined, using each of the above models, by a method dependent only 
upon the physical characteristics of the chosen atmospheric model. Placing 
now additional weight on the results obtained using the Swihart model, the 
most probable value of the nitrogen abundance is found to be log Ay = 7-86 
(on the basis log Aq = 12-00). 





Introduction.—In an earlier paper (9), concerned with the determination of 
solar atmospheric abundances, a method of determination was developed, 
the Planckian gradient method, which was based entirely upon the physical 
characteristics of the assumed theoretical model atmosphere. In order to examine 
the numerical predictions of the method, and particularly to illustrate the diver- 
gencies when compared with the method of weighting functions, an element of 
high ionization potential was considered (nitrogen, I.P. 14-54 volts). It was 
also of interest to examine the sensitivity of the method to the particular atmos- 
pheric model and therefore three models were chosen: (i) the model adopted by 
Claas (2), which is referred to as the ‘‘ Claas model ’’, (ii) the model of Vitense 
(14) and (iii) the Swihart model (12). The method of derivation, and the 
detailed models themselves, are described in the aforementioned paper. It 
was found that the derived nitrogen abundance was moderately insensitive to 
the adopted atmospheric model. In obtaining an estimate of the most probable 
value of the nitrogen abundance, however, it is necessary to obtain some measure 
of the reliability of each of the atmospheric models. It is the purpose of the 
present paper to obtain such a measure. 

The criterion of reliability adopted here is the correct prediction, by the 
model atmosphere, of the observed line profile in the far wing of the Hy absorp- 
tion line (A= 4340°5 A). The Balmer lines are chosen since they are very strong 
and can be measured accurately in the far wing; in addition the excitation potential 


14 
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of the lower level (10-20 volts) is very close to the excitation potential of the nitrogen 
lines (10°33 volts). A consequence of this latter point is that the Balmer lines 
have approximately the same temperature sensitivity as the nitrogen lines and 
therefore originate in the same regions of the model atmosphere. 

Of the various physical mechanisms producing a broadening of spectral lines, 
e.g., natural broadening, Doppler effect, collision damping, Stark effect, etc., 
it is found that, in the far wing of strong lines, the contribution to the absorption 
produced by the former effects is negligible and line broadening is almost entirely 
due to the Stark effect. We wish then to predict the line profile in the far wing 
of Hy, i.e., sufficiently far from the line centre that effects other than the Stark 
effect may be neglected but not so very far in the wing that observational measure- 
ments become uncertain. 

1. Kolb’s theory of the first order Stark effect.—Most theoretical considerations 
of the problem of Stark broadening, outside the narrow Doppler core, have 
been based upon the Holtsmark (4) statistical theory where the radiating atom 
is subjected to long-range coulomb forces due to neighbouring ions which produce 
an ‘‘average’’ field. The slowly moving ions are treated as though stationary 
during the times of interest and the field produced is then dependent only upon 
the ion density. Kolb (5), however, has refined the theory and takes into account 
not only the ions but also the effect of the electrons and their finite velocities. 
A detailed account of Kolb’s theory has recently been published by Aller and 
Jugaku (1); in the present paper we merely restate Kolb’s results for the far 
wing of the Hy line while retaining Kolb’s original notation. 

|! In Kolb’s theory the ion broadening is treated according to the usual Holts- 

mark theory, where the ions produce an average field at the radiating hydrogen 
atom which produces a Stark resolution of the degenerate hydrogen substates. 
The high velocity electrons broaden these static field states since (a) they produce 
phase shifts in the emitted wave train according to the ‘‘ weak collision ’’ theory, 
and (b) they induce transitions between the degenerate states (the non-adiabatic 
effect). The total absorption coefficient will be proportional to the sum of the 
individual broadened Stark components weighted in accordance with their 
theoretical intensity. In the present theory for Hy the non-adiabatic effects 
have been taken into account in only a very approximate way, i.e., Kolb has 
calculated the non-adiabatic effects accurately only for Lyman-« (where he 
found that the adiabatic contribution was approximately half that of the non- 
adiabatic effects) and he has tentatively assumed that approximately the same 
relative contribution holds for the Balmer lines. 


If Aw is the displacement, in circular frequency units, of the Stark component 
from the line centre, i.e., 


Aw=2nAv=2nc (1) 


(for Hy, A= 4340 A, Aw=AA.10™) we may define a parameter 


Aw 
, (2) 


as a measure of the ion field, where 7, is the ion half-width for the Holtsmark 
average field strength: 


Vu=4°524n,X;,. (3) 
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Here n, is the ion density (which is taken to be equal to the electron density n,) 
and X,, is an integer which is characteristic of the displacement of the kth Stark 
component. Similarly we define a parameter 


Aw 


7 << (4) 
Yk 


as a measure of the electron field, where y, is the half-width for electron broadening. 
Kolb finds 


n, 

Y¥,= 122 x 10 +X} 
VT 

where F, is a function of the parameter 5, and 


/ 
5, = 1146 x 10-*X, +, (6) 


actually 
F,, =0°2124 — 0°383 log5,. (7) 


A factor 2 has been included in equation (6) to allow, in an approximate way, for 
the non-adiabatic effects which were calculated accurately only for the Lyman 
series. 
The line absorption coefficient per absorbing atom for each Stark component 
k is defined as 
14(Ato) d(Aw), (8) 


with the normalization condition 


y an x,(Aw) d(Atw) = an( = =) fo (9) 


where f,, is the total oscillator strength of the displaced components (for Hy, 
f, =0°0377). Weighting each component by its theoretical intensity Kolb finds 


1y( Aw) d(Ato)= an (= =) fo 2 fs ~ I,(Aw) d(Aw), (10) 


where f, is proportional to the f-value of the kth component, and /,(Aw), the 
intensity of the individual Stark components, is given by 





x8 5A 
a ee 


Here the parameter A,(x) has the value 





_ 1 sin(7/2tan-' x) 
af I Saori (1 + 1/x?)"/4 (12) 


which has been tabulated by Kolb (5). In order to obtain the total absorption 
coefficient at any value of Aw we must now sum over all Stark components: 


2(Aw)= Say(Ae)=[ a (7) fo? |Z $5 In(Au), (13) 
14* 
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'‘ Introducing the numerical values of the constants and writing s,=f,/>f, we 
have (noting equations (2) and (11)) 


a(Aw) =0-00400 S's, — I, (Aw) 
k Yk 





=o: x x? 1:175A,(x) 

o-00400 Sau | es + (Bs) et +] 
oe 1:175A,(x) (7, \*? 

= 24 at Ansa (5%) X,37+.. | 


17 2 
“— oo 254s (x)+ 2 wie (%) r5eX? 24i(x)+.... (14) 
" 


It is found that the higher terms are unimportant in the far wing, i.e., 
when Aw is large compared with the half-width for electron broadening 


(= do, :) When this condition holds 
Yk 





Mihi date > (15) 


(x? +1 tif x 
and the equation, valid for the far wing, becomes finally 


' I , 5, \3i2 em 

a(Aw) = 2 oe 25+ eke (%) S.X,?*As(x). (16) 
The second term in equation (16) is the usual Holtsmark expression (see, e.g., 

; Aller and Jugaku (1)) with the exception of the factor 4,(x) which, however, is of 
the order of unity for large values of x. It will also be noted that the absorption 
' coefficient is dependent upon both density and temperature, unlike the result 
obtained by the Holtsmark theory which is temperature independent. 

The Kolb theory, essentially in the above form, has been applied to A star 
atmospheres by Osawa (ro) and to B star atmospheres by Elste, Jugaku and 
Aller (3); these investigators conclude that the Kolb theory predicts the observed 
line profiles more closely than does the Holtsmark theory. More recently the 
theory has been revised and generalized by Kolb and Griem (6) in a form applicable 
to the line centre as well as to the wing. Kolb (7) points out that the revised 
theory indicates that the shape of the distribution is the same as that given by the 
above theory, although the required correction to the Holtsmark theory may be 
slightly larger than indicated above. 


2. Calculation of the line absorption coefficient 1,—According to Kolb the 
f-numbers for the Stark components of Hy are 


Xr fe Xe Sk Xe fk 

° 1416 10 1760 20 80 

2 156 12 166 22 7 

3 922 13 1664 

5 192 15-1152 DLS = 9016 
116 17 52 

8 15 18 1318 


the factor s, may then be calculated directly. 
The calculations are carried out in the following steps: The atmospheric 
model provides details of the electron pressure P, (and hence n, since P,=n,kT) 
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and temperature as a function of optical depth, i.c., r=7(P.,T). For given 
values of Aw and + we may compute: 


(i) 5, from equation (6) and F, from equation (7) for each Stark component k ; 
(ii) X,2F;, and hence y, from equation (5); 
(iii) 7,/X, from equation (3); 
(iv) the factors, variable with k, of the 1st and 2nd terms of equation 
(16) for each Stark component k; 
(v) the sum of (iv) over all Stark components ; 
(vi) « in accordance with equation (16). 


This process was repeated at moderately regular intervals of optical depth, 
and at three distances from the line centre AA=5, 6, 7 A corresponding with 
Aw=5, 6, 7, x 10! (see equation (1)), for each of the three model atmospheres. 
It is found that the ionic and electronic contributions to the absorption coefficient 
are of the same order of magnitude and therefore, in the far wing, the electronic 
component should not be neglected. 

The line absorption coefficient per absorbing atom, «,, has thus been obtained 
in accordance with the Kolb theory and the transition to the mass absorption 
coefficient /, may now be made by the relationship 


N; 
l= Numa his (17) 


where N; is the number of absorbing atoms per cm*, and N,m, gm cm-~* is the 
assumed density of the solar atmosphere (due to the overwhelming abundance 
of hydrogen). Introducing Boltzmann’s equation we have 


N,; .. &.10-"° 5040 
ek 6b Wileeeenss 6= 2 8 
Nu u ‘ T (x ) 


where g; is the statistical weight (2n*) of the lower level concerned in the transition 
and V is the excitation potential of this level (V = 10-196 volts). The partition 
function for hydrogen has been calculated by Claas (2), i.e., u=2-0. Equation 
(18) then becomes 
: en 
L, es my Hy (19) 

The line absorption coefficient /, can therefore be obtained for the optical depths 
7,=7,(@) at which «, has been previously calculated and graphical interpolation 
provides /, at regular intervals of r,. 

3. Prediction of the wing of the Hy absorption line.—The intensity dip in 
the continuum, at the centre of the solar disk, for radiation of frequency v is 


°dB, Tvl, 
i { “ae [ 8 it, |ar, (20) 


for a line formed by pure absorption, where B,(@) is the Planck function and 
k, is the continuous absorption coefficient. Equation (20), used by the author 
for abundance determinations in a previous paper (g), is found to be unsuitable 
for the present purpose since the ratio /,/k, is very small, except at large optical 


depths, and the integration | l,/k, dr, is therefore rather uncertain. An alter- 
0 
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Fic. 1.—Contribution to absorption as a function of optical depth for the wing of H-gamma 
(5, 6 and 7A from line centre). 


native form of this equation (see Mugglestone (9) equations (6.4) and (6.5) 
there) previously obtained by Unséld (13) is preferred, i.e., 


ar,={"| ” Byewrdr, Bye" | dr, (21) 
0 = k, 
which was numerically evaluated in the following manner :— 


(i) P(r))= [ B,e-re Edy B,e-*». 
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The ratio of absorption coefficients k, (A= 4340 A) to ky (absorption coefficient 
at standard wave-length A= 5000 A) were obtained from tables supplied by Dr B. 
Strémgren (unpublished), hence the relationship between the corresponding 
optical depths could be calculated from 


T,= ["jeare (22) 


0 0 
at regular intervals of optical depth rt». The particular atmospheric model adopted 
provides t,=7,(@), and since for a given wave-length (A= 4340A) the Planck 
function is a function of optical depth only, i.e. B,=B,(@), B,e~"*k,/ky may be 
determined for each value of 7). The required integrations 


” k " k 
[. dite dry— {" le 


were calculated in accordance with the Gaussian quadrature formulae and hence 
the function P(r.) is obtained at values of 7) with regular interval. (See Fig. 1.) 


(ii) AT = [EPC dr). Values of the continuous absorption coefficient 
0 0 


ky=k,(0,P.) are available from tables of hk /k=k,/k(0,P.) (Strémgren, 
unpublished) and k=k(6,P,) (Strémgren (11)), and the line absorption 
coefficient /, is calculated directly from equation (19) utilising the Stark effect 
results obtained from Kolb’s theory. The results obtained for the integrand 
L,/ky.P(79) are shown graphically in Fig. 1 and the resulting AJ,, calculated 
in accordance with the Gaussian quadrature formula and expressed as a percentage 


of the intensity in the continuum, are given in Table I, i.e. the intensity in the 
continuum J¢, which is given by 


Is= [” Byer (23) 
0 


0 
has been computed in section (i) and hence the predicted AJ,/J° may be compared 
with the observed intensities in the wing of the Hy line as measured in the Utrecht 
Photometric Atlas (8) (see also Fig. (2)). 


TABLE I 


Predicted dip in the continuum AI, expressed as a percentage of the continuum intensity 
(i.e. AI,/I®) and compared with the observed intensities in the wing of the H-gamma line. 
Model 

atmosphere AA=5A AA=6A AA=7A 

per cent per cent per cent 
Claas model 14°1 9°3 
Vitense model 13°6 8-9 
Swihart model 10°5 6°8 
Observed 8-8 6°3 


4. Conclusions.—The results obtained in the above calculations, together 
with the results of the recently revised theory of Kolb and Griem (i.e., that the 
electronic contribution to the absorption coefficient should be further increased, 
thereby increasing the predicted absorption in the far-wing), clearly indicate 
that, on the basis of this criterion of reliability, the Swihart model should be 
preferred. It is clear however that this preference should be applied only for 
absorption lines formed mainly in the lower layers of the atmosphere and for 
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Fic. 2.—Wing of H-gamma. Intensities predicted by theoretical models compared with observed 
values. 


which equations (20) and (21) are applicable: i.e. weak lines (or the far wings 
of strong lines) of high excitation potential. 

With regard to the nitrogen abundance in the solar atmosphere a previous 
paper (9) established the following results: log A,=7-71 (Claas model), 
log A, = 7°81 (Vitense model) and log A, = 7-93 (Swihart model). In view of 
the results of the present paper it has been decided, somewhat arbitrarily, to 
assign triple weight to the results obtained using the Swihart model and single 
weight to the Vitense and Claas models. On this basis the most probable value 
of the nitrogen abundance is 

log A, =7°86 (log A, = 12-00). 

The author would like to thank Dr A. C. Kolb of the United States Naval 

Research Laboratory, Washington, for advance information with respect to his 


theory of Stark broadening, and for private correspondence entered into in 
this regard. 


Department of Physics, 
University of Queensland, 
Brisbane, 
Australia : 
1959 August. 
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TRANSITION PROBABILITIES FOR FORBIDDEN LINES OF Ne Iv 


R. H. Garstang 
(Received 1959 September 18) 


Summary 


Transition probabilities have been computed for magnetic dipole and 
electric quadrupole radiation for transitions between the levels of the 2p* 
configuration of Ne tv. The results of a calculation of the relative intensities 
of lines of the *P-*D multiplet are compared with observation. 





The forbidden lines of neon are of particular interest in planetary nebulae, 
not least because lines of three consecutive stages of ionization (Ne 11, Ne Iv and 
Ne v) are frequently observed in the same object. Published lists of transition 
probabilities usually include a few lines of [Ne 1v], but for calculations on the 
physical state of the nebulae information is needed on all transitions in [Ne Iv], 
whether in the accessible spectrum or not. ‘The purpose of this paper is to 
provide such a complete list. ‘The calculations involved are similar to those for 
[O 1], described by Aller, Ufford and Van Vleck (2) and by Seaton and Oster- 
brock (12). We use the same notation as in their papers. 

The spin-orbit parameter {.—'This was estimated by interpolating (1/4 between 
Ne and Nev, Natv and F 1v, and Mg v and O11, the values of ¢ used for all 
these ions being those of Garstang (§). The values of { obtained by the three 
methods agree very closely, and the fina] value adopted for Ne Iv is =692 cm™?. 

The spin-spin parameter 4.—-Interpolation of 7! between Ne m1 and Ne v, 
Na Iv and F tv, and Mg v and O 111 gives, respectively, 7 = 4°92, 4°28, 4°17 cm™}. 
Examination of the values of 7 for all these ions (§) shows that the value of »!* 
for Ne 111 lies off the best linear law of »' against Z for the 2p* configuration. 
If the value of 7! is smoothed before interpolation between Nem and 
Ne v a value »=4:49 cm ! is obtained for Netv. As a further check » was 
computed from the self-consistent field with exchange wave function of Froese 
and Hartree (4), the result being »=4-48cm''. ‘Taking all these values into 
consideration 7 = 4:40 cm™! was adopted for Ne Iv. 

The term separations (PS), (PD), (DS).—Using energy levels given by 
Moore (8, with additions in Appendix to Volume III based on Bowen (3)), we 
adopt (PS)=62117 cm™!, (PD)=21140 cm ', (DS) = 40977 cm". 

The quadrupole integral s,.— This was estimated by interpolating between the 
values of s,-!/* for the 2p? and 2p* configurations (§). This gave 

$= i r? P*(2p) dr = 0-273. 


0 


As a check, the wave function of Froese and Hartree (4) was used, giving 
5, = 0273; this value was adopted. As a further check, the value of s, for Ne v 
estimated (§) from other ions in the 2p* series is 0-235; Froese and Hartree gave a 
wave function for Nev and using it we obtain s,=0-240. The excellent 
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agreement leads us to think that our estimate of s, for Ne Iv is as satisfactory as 
can be obtained at the present time. 

Doublet intervals—These were computed using formulae (16) of Aller, 
Ufford and Van Vleck (2). The results are listed in Table I, together with 
those of Horie(7) and the most recent experimental values (8, Appendix to Volume 
III). The agreement we obtain is quite good, and is rather better than that of 
Horie. In order to examine whether the results are sensitive to the adopted 
value of », the calculations were repeated using 7 =4:80 cm", all the other 
parameters remaining unaltered. ‘The doublet intervals are given in Table I. 


Vol. 120 


TABLE I 
Doublet intervals in ?P and *D terms 
AP) AD) 
(cm~') (cm~') 
This paper (7 = 4°40) +3°0 — §3°1 
Horie (7) —3 —65 
This paper (7 = 4°80) —o'8 


—61'°5 
Experimental (8) +64 


— 44°'8 
It is seen that the agreement between calculated and experimental values is much 
better for y = 4-40 cm™' than fory = 4:80 cm™!. Indeed, 7 = 4:00 cm ! would give 
nearly perfect agreement for the doublet intervals. However, because such a 
value would 'be well below that obtained by interpolation and from wave func- 
tions, we prefer to use 7 = 4°40 cm™ in further calculations. 

Transition probabilities.—The transformation coefficients expressing the 
perturbed wave functions in terms of the unperturbed wave functions were 


obtained from equations (18) of (2) supplemented (12) by a” =1, 6” ={/(PS) and 


ce” =(4/5)¢/2(PD). The transition probabilities were then obtained from 
formulae given by Shortley, Aller, Baker and Menzel (13, equations (2), (3), (11) 
and (12)). Our final results are given in Table II. Comparison with the few 


Tasie II 
Transition probabilities for [Ne tv} 


Lower 
*P iia 
‘Djs 
"D2 
D5 


*Ss2 


*D,: 


1°33 
0°53 
1°4 x 107% 


A, 
26 x 107% 
o'1g 
7°9x 107? 
oll 
o'16 
15 xX 1077 
8-6 x 1078 
1'rx107!7 


A total 
2°3 x 107° 
0°40 
0°44 
oll 
0°39 
1°33 
0°53 


1'4 x 107% 
59x 107" 
56 x 107% 


‘Sais 18x 1074 


4°1 x 1074 
*Das Sas 5°3 x 107° 


2:7 x 10-* 
transition probabilities published by Pasternack (11) shows differences mainly due 
to revised doublet intervals and our lower value of s,. Comparison with Naqvi’s 
results (9) shows differences which may be ascribed to our use of empirical 
term separations; reasons have been given (6) which lead the writer to the view 
that it is better to use empirical term separations instead of attempting to deter- 
mine the Slater integral F,, as was done by Naqvi (cf. Naqvi and Talwar (r10)). 
The calculations were examined to see whether any of the transition pro- 
babilities were sensitive to the adopted values of the parameters. The quadrupole 
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parts of the ?P;.—*Ss). and *P, ,.—*S, . transitions appear to be somewhat sensitive 
to the value of 7, but as the magnetic dipole components predominate for these 
two lines, the sensitivity has no practical consequence. None of the other 
transitions depends critically on the adopted parameters. 

Intensities in the*P—*D multiplet.—The relative intensities of the components 
of this multiplet have been estimated by Aller, Bowen and Minkowski (1) in the 
spectrum of NGC 7027. Using the transition probabilities given in Table II 


Tasie III 
Observed and calculated intensities in ?P —*D multiplet 
Wavelength Transition Observed Calculated 
4714°26 *P52—"D5,2 o'8 vr 
4715°65 ‘ 12> *Ds)2 o's o"3 
472416 *P32—*Ds5)2 1°3 13 
4725'60 *Pi 2 —* Daj 13 12 


and his own estimates of the collision strengths, Dr M. J. Seaton has calculated 
the expected intensities of the lines of the ?P-*D multiplet. His results are 
listed in Table III; the theoretical intensities are normalized to the same total 
as the observed intensities. ‘The agreement is well within the errors of calcu- 
lation and observation. 

Acknowledgment.—'Thanks are due to Dr Seaton for suggesting this work, 
for contributing Table III, and for his continued interest in the writer’s calcu- 
lations. 


University of London Observatory, 


Mill Hill Park, London, N.W.7: 
1959 September 14. 
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CAN SPHERICAL CLUSTERS ROTATE? 


D. Lynden-Bell 


(Communicated by R. v. d. R. Woolley) 
(Received 1959 September 8)* 


Summary 


It is shown that contrary to former belief an unrelaxed spherical cluster of 
mass points may rotate without becoming oblate. When used of an 
unrelaxed system the argument that spherical objects do not rotate is thus false. 
Some present estimates of the Sun’s circular velocity round the galaxy rest 
on this argument applied to the system of globular clusters. 





Introduction.—We are here discussing clusters of gravitating mass points. 
Thus applications may be considered at several levels. 


(i) Clusters of galaxies. 
(ii) Fo galaxies themselves. 


(iii) The system of globular clusters centred on the galactic centre. 
(iv) The globular clusters themselves. 


The system of globular clusters shows no pronounced flattening at the 
galactic poles; is this good evidence that it does not rotate? If it is not, then 
the justification for the determination of the Sun’s circular velocity from the 
Doppler motions of the globular clusters is destroyed. These observations give 
the Sun’s velocity in the coordinates in which the system of globular clusters 
does not rotate, and we have no longer any justification for assuming that this is 
not a rotating coordinate system. 

I wish to show that unrelaxed systems may be steady, spherical, and rotating. 
Perhaps there are cosmogonical reasons why a spherical system may not rotate, 
but such a conclusion cannot be reached on dynamical grounds alone. It is 
possible that some systems with net angular momenta evolve into spherical 
states in times short compared with their relaxation time, before eventually 
flattening to the relaxed ‘‘gaseous’’ forms discussed by Jeans (1,2). The 
form of the frequency function for the general spherical cluster is found in this 
paper. In general these clusters rotate. This is in contradiction to Jeans’ result, 
but is obtained by using his method correctly and following the consequences. 

My deepest thanks are due to Dr Woolley who not only suggested the problem, 
but convinced me of the result. The argument that follows arose from stimulating 
and amusing discussions with him. 

Consider a mighty Maxwell demon whose job it is to reverse the directions 
of motion ef stars so that they describe their orbits in the reverse sense. We shall 
assume him capable of violating the law of angular momentum conservation. 


* Received in original form 1959 May 6. 
+ The problem came originally from Dr Sandage. 
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On his entering a non-rotating spherical star cluster let him pick a direction and 
a sense of rotation about it. Let him reverse the velocities of all those stars which 
do not have the picked sense of motion about this axis. It seems likely that he 
can do this without disturbing the formerly spherical density distribution. If this 
is the case he succeeds in making a spherical but rotating cluster. The assumption 
implicitly made is that any non-rotating spherical cluster may be considered as 
two inter-penetrating clusters, rotating in opposite senses, each of which is itself 
spherical. (it is clear that they cannot be flattened since then the complete 
cluster would also be flattened.) 

Jeans’ theorem (1, 2, 3, 4).—This states that for a steady system of mass 
points moving independently in a potential V, the frequency function* f(which 
gives the density of representative points in the 6-dimensional phase space) 
must be a function of the independent first integrals of the equations of motion 
only. When V is spherically symmetrical 4 such integrals, the energy and the 
angular momenta, are readily obtained. A fifth may also be obtained but it may 
be shown that continuity of the distribution function forbids dependence on this 
integralt. We shall therefore neglect it. 

We write for the position, velocity, velocity moment, and energy: 


r= (x,y,z) Cartesian = [r, 0,4] spherical polar. 
c = (u, v, w) Cartesian = {c,0,,,}. 9,.,¢, give the direction of c. 


W=6 x C=(w,,a2,H3)=[w,4,,4,,]. 9 bey give the direction of w. 
e=c/2-—V. 


We also use the notation } ... dc for | | | ...dudv dw and | ... dr similarly. 


Jeans’ theorem reads 
f=f(@, ). (1) 


We assume that stars of different masses are distributed with similar frequency 


functions. This gives the same result as assuming all stars to have the same 
mass m1. 


The potential V may be made up of two parts; 
(1) Vi the potential due to the cluster of moving mass points; and 


(2) V4 corresponding to an external field applied to the system. 
Poisson’s equation reads: 


V°V seis = — 47 YP seis = — 47M [fee (2) 


Spherical clusters 


Jeans argues that a spherical cluster is invariant under rotations. Hence f 
must be invariant. The independent invariants of w and ¢ are w® and ¢ and 
hencet 


f=f(w, «). (3) 


* The cluster is completely described by its distribution in phase space. The density of this 
distribution is given by the frequency function which is often called the distribution function. 

+ A paper is in preparation on the justification for the neglect of the further integrals, such as 
this, in the use of Jeans’ theorem (3). 

t By the theorem that an invariant function is a function of the invariants. Our result may be 
obtained elementarily by partial differentiation (5). 
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However the condition that a cluster be spherical is not that f should be 
rotationally invariant but that » should be. The problem of finding the most 
general frequency function for spherical clusters thus entails finding the most 


general function f(w,«) for which £ = | f d®c is a function of r alone. 
m 


There is one further condition; f must be nowhere negative. We shall 
call a cluster of form (3) a Jeans cluster. If also f is nowhere negative we say it 
gives a ‘‘real’’ Jeans cluster and if conversely an ‘‘ imaginary ’’ Jeans cluster. 

Our first problem is that of fitting the density of a Jeans cluster to that of a 
given spherical cluster. 

Theorem.—To any given spherical cluster (with frequency f,) there corre- 
sponds a Jeans cluster (f,) of the same density. 


Proof.—Define 


fr =f; = — i fi(@, «) sin 6, d6, dd. 
47 
Clearly f, is a function of (w?,e) only and it thus gives a Jeans cluster. It is 
also non-negative (since f, was) and therefore gives a real Jeans cluster. It is 
obtained from f, by averaging over the directions (#,,¢_,) of the normals to the 
planes of the stellar orbits. This averaging process merely tilts the orbital planes 
so that these normal directions are evenly distributed. This process does not 
change the’ mass or the speed distribution* in the spherical shell r to r+ dr. 
Thus the radial variation of density is the same in f, and f,. Since they both give 
spherical clusters the density only varies radially and they give the same density 
everywhere. 
We define the deviation Af by: 
Mf=fhi-fe 


By our theorem 
| Af d¥cxo all r. (4) 


The general Jeans cluster is known, thus the general spherical cluster may be 
found if we can find the most general form of the deviation Af. 
We wish therefore to find the most general Af satisfying (4) and subject to 
the restriction 
fr,+Af2o. (5) 
This last may be considered more as a restriction on the type of Jeans cluster (f,) 
which may be added to a given deviation than a restriction on the deviation itself. 
We therefore consider deviations satisfying (4). 
A particular class of solutions.—Consider a Af antisymmetric in @, i.e. 
Af(es, ¢)= — Af(—a, ¢) 
@=rxc, € is even in c, thus for any given r this Af is antisymmetric in ec. 
Hence (4) is satisfied. If we add an f,(m*,«) such that (5) is satisfied we will 
have succeeded in making a non-Jeans spherical cluster. Note that creation of 
such a Af antisymmetric in w is just what is achieved by the demon on reversing 
stars in their orbits. In general he has not been given freedom to turn all the 
* Speed =| velocity |. The result will be needed later. 
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stars one way, but is no doubt chewing his claws with discontent over not being 
allowed to finish the job. Later we shall find clusters which cannot be obtained 
by letting a demon act on a Jeans cluster. 

We take as a model the Jeans cluster due to Schuster (6). 


(1) Frequency function 
7/2 
fa=A (v- <) 


1.48 @M/( 1) 
“74/273 m\yM) * 


(2) ‘Total mass M, mean radius 7 = 2a. 
(3) Potential 


where 


Density 


(4) R.m.s. stellar velocity 


Je S V35( 2)". 
4 2a 


(5) Average magnitude of the angular momentum of a star taken about any 


mam» (2 a) / 
7 2a 


axis 


Total angular momentum of cluster H =o. 
Consider the deviation given by 
lh 
Aifice niga Se 
f 1+ (w,/h)* ’ 
where h>o is a constant. Since the minimum of 5 } hence 
I+x 
fo+Afzo. Thus, since Af is antisymmetric in w, f;=f.+Af is a cluster with 
the same density as f, everywhere. Its angular momentum is 
° ” led a. 
H= | [rah de ar | [fers ae dtr =n | eps Pedr>o. 
I 
We have succeeded in making a rotating cluster. A different cluster for which 
the integrals are tractable is the demon’s cluster. ‘This however has a 
discontinuous deviation 
=+I, @M,>0, 
Af=Sgn(a;)f. where Sgn(w;) 
=-I, Ws, <0. 


The properties (1)...(4) hold for this cluster too but the total angular 
momentum 


H= i | [ler faded 


since w, Sgn w;=|m,|. On integration 


Baa) \ a) = a) ae 
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On comparison with (5) we see that the contributions from the individual stars 
all add, rather than subtract, as the demon requires. 

The general spherical star cluster.—We shall need the following mathematical 
theorem. 

Theorem.—The most general function defined on a sphere, whose integral 
around every great circle vanishes, is any function antisymmetric on inversion 
in the sphere’s centre, i.e. if the sphere is |?]=1 


a(P)= —g(—P). 

To prove this theorem we need the concept of the great circle transform 
which I will now develop. 

To each pair of antipodes r, —r on the sphere there corresponds one great 
circle, their equator. ‘Thus the space whose elements are great circles is similar 
to the space of antipodal points. ‘lo make up for this pairing of the points we 
shall until further notice only deal with functions on the sphere which are 


symmetric under inversion. ‘The great circle transform of a function f(@,¢) is 
defined by 


Se ’ 
He B)= = |, Jb) dss 


where ds, is the infinitesimal path length round the great circle C,, which has 
(a, 8) as: one of its poles. Since the choice of pole is arbitrary f is symmetric. 
Thus we have a transformation on the sphere which, to any symmetric function, 
yields a corresponding symmetric function. We shall here assume that there 
exists an inverse transformation* such that given f we can find f. 

Suppose f is the symmetric part of a function whose integral round every 
great circle vanishes. ‘Then from the definition of the great circle transform 
f=o. However since there is an inverse transform this defines f. In particular 
since the transformation is linear we have f=o0. Clearly the antisymmetric part 
of the solution, g say, can be anything since a great circle through P will always 
pass through P’s antipodes. Hence the theorem as stated. 

We now use our theorem to prove the result for star clusters. We left this 
problem at equation (4) where we reduced it to the problem. What is the most 
general function Af(@, «) satisfying 


for d*c=0. 


Measure 6,, ¢,, from f as pole c=(c,6,,¢,) we require 


| | fay (re sin 0, @, & - V(r) sin 0, d8,c? de dg, =0 


where oe ee 

hed 
note & is independent of |r|, |c|, @.. It merely gives the direction of the normal 
to the r, c plane (in 3 space). We may thus perform the first two integrations, 
defining 


ee 


2 
g(7,@)= | far(resin 8.5, — _ v) sin 6, d0..c? de 


and obtain 


[etre a8, =o allr. 


* I have a rigorous proof of this result but it is neither short nor beautiful. 
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Now @ is the normal to the plane ¢,=const whereas all such planes pass 
through #; thus this last integration is in effect over all directions & perpendicular 
to a given #. Since r is not involved in the integration, regard it as a fixed 
parameter. Then translating (6) into words :— 

A function of direction is such that when integrated over all directions 
perpendicular to a given # the result vanishes. However the given # is arbitrary; 
thus g(r, @) is a function of direction such that its integral over all directions L ” 
to any given direction vanishes. Alternatively g is a function defined on a sphere, 
whose integral round any great circle vanishes. Hence by our theorem 


&(r,&) = —g(r, —&) (7) 
i.e. the symmetric part of g vanishes. 
Split Af into symmetric and antisymmetric parts writing Af=A,f+A_f 
where 2A, f=Af(w, «)+Af(—@,«), (7) implies 


{| A, fsin@, cdc= o. 


Since & is not involved in the integration we may now regard it as a fixed 
parameter. Write A,f=A,f(@;m*,«). (7) declares that, for each &, A,f is 
the difference between two Jeans distributions of the same density. There are 
many Jeans distributions (f(w*,«)) corresponding to a given density function* 


or) =a { f@c= ant | | fsin 0, d0.c? de. 


Any distribution of the form f,+A.,f may thus be obtained by having different 
corresponding Jeans distributions for each coplanar set of motions (®.) 
Finally we obtain the general distribution by adding any antisymmetric A_f 
provided the total f,=f,+4,f+A_f is nowhere negative. Summarizing our 
method we have the following steps. 


(1) Choose a density function p(r); Poissons equation (2) yields V. 
(2) Find the Jeans distributions corresponding to that density*. 


(3) Consider a function f,(@,w?,«) which for each & reduces to one of the 
above set of Jeans distributions. f, may be found by assigning one of 
this set to each direction &. f, gives the general spherical cluster 
distribution symmetric in @. 

(4) fi=f,+4_f is the general spherical cluster distribution, where A_f is 
antisymmetric in w and is restricted to make f,+A_f20o everywhere. 

In short, take the Jeans distributions of the right density; assign one such to 
each orbital plane. Let the demon reverse some motions, and we have a way 
of constructing the general spherical cluster. 

The fastest rotating clusters 

‘* Given the density p(r) of a cluster what is the greatest angular momentum 
it can hold consistent with its maintaining that density ?’’ 

‘* What is the spherical cluster of greatest angular momentum irrespective 
of density distribution ?”’ 

This latter is rather a wonderful figment of the imagination; a great spherical 
shell of stars round which each star crawls with a circular velocity which is quite 
minute. However all their senses of rotation are the same about the demon’s 

* See Appendix. 


15 
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axis, and though their speeds are small they carry large angular momenta. We 
now return to the saner question posed above. We first demonstrate the 
existence of spherical clusters in which all the stars move in circles. By 
considering the orbits in each plane we then show that circular orbits can carry 
more angular momentum about a given axis that any other distribution of orbits 
that give rise to the same density. This is sufficient to demonstrate our result. 

We shall assume the given density distribution to be such that the circular 
orbits are everywhere stable. 

Since p(r) is given, hence V(r) is known, hence r(e) (the radius corresponding 
to circular orbits of energy «) is known; and w[e] (the angular momentum in 
circular orbit of energy €) is known. 

It is a well-known property of the stable circular orbit that ‘‘ Of all the orbits 
of given energy the circular orbit carries the greatest angular momentum w’’. 

Construct a Jeans cluster as follows: 

f(w*,e)=0 unless w*?=(a[e])’. 
This ensures that all the stars move in circular orbits. By giving f((w[e])’, €) 
as a function of « we may populate each spherical shell just as we wish. In 
particular we do this so that we obtain the required density 
p(r). (8) 
We have thus proved that: 

‘“*'To each cluster there corresponds a Jeans cluster with the same density 
everywhere and in which all the stars move in circles.”’ 

Consider the planar motion of a star. This motion gives rise to a certain 
distribution of density with radius ‘‘p(r)’’ and a certain angular momentum 
mh cos \ about an axis at A to the plane’s normal. Can we replace this star by a 
distribution over orbits which taken together give the same density p(r) but a 


greater angular momentum ewcosA? In particular does replacement with 
circular orbits have the desired effect ? 


_— : . 
5t= — gives the time spent by the star in r>r+dr. Let 7,,, and 7, be 
; 


the limits of this star’s orbit. Thus , 
ar 


r 


(‘Tmax dr 
|, 


io F 
« "min 
gives the ‘‘number of stars’’ which, permanently within r-r+dr, will give 
the same average density there. Their angular momenta is given by 


tiver(r) va cos A 
? 


"dr 
17 


where #im(r) is the circular angular momentum at distance r._ The total angular 
momentum gained by such replacements is 


ee am(r) on cos A 
? 


¢ ae —hcosa 
maz dr 





'min * 
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The replaced system thus carries more angular momentum than the original 
provided (calling the integral J) 


Joos= |" (w(r)— h) = cosd>o, 
but min 


72 
oy ei Vek 
2 

and 


(w(r))?__ aV 


rs or 
thus the integral J is 


a 


Tmax, po dr 


/(—V)—h 
| ine (v )Taet V)—hejr2 


= (— V’ —h?®/r°) dr 
min 7-32(/ — V' + h]r®?)V 2(E + V)— her 
Integrating by parts (V 2(E+ V)—A?/r? vanishes at end points) 


(Cr. 
J=o+ | = 


min 




















V2EFV) =e (WV + h/r5h)-4732) dr, 


therefore 








[ maz \/2(E+V)—hr 


7 (Jr thir) pl a =)- sr |e 


#2) vey M 
dr re) dr eS 


d . 
= (-V'r)<o 


(circular orbit stability) and V’<o always. Thus 


¢ Ie _-# <o 
dr =) e : 


hence J>o. We may take A<7z/2 by choosing the sense in which the stars 
describe their orbits. Thus JcosA>o. We have therefore shown that angular 
momentum is always gained in the replacement of non-circular orbits by circular 
orbits giving the same density, i.e. the orbits carrying greatest angular momentum 
for given density are circular and the net angular momentum is greatest if the 
senses in which they are described are arranged by the demon. Thus: 

‘““The fastest rotating spherical cluster of given density is the ‘circular 
demon cluster ’.’’ 

It is interesting to note that if all the stars travelled at the same speeds, at 
the same central distances, and all contributed their angular momentum in the 
same direction instead of at all angles up to 7/2*, then the total angular momentum 
would be twice that of the above cluster. This indicates how little is the relief 
granted to a rotating cluster when it flattens. 


15* 
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APPENDIX 
How to find all the Jeans distributions of a given density: 
Find one such distribution, e.g. construct the Jeans cluster in which all the 


stars move in circles as in (8)*. Let its distribution of speeds with radius be 
F,(c,r), i.e. 


ign | Fa(c,r) de=p(r). 
Take an arbitrary function A’ F(a, r) satisfying 

[aru r)demo alls 
write g=c and A’ F(e,r) =AF(c,r) then 


[are dc= o. 


‘{ 
F, = F,+ AF has the same density as F, and all the F with this property may be 
so obtained. We show that to each such F, there corresponds a Jeans distribution. 
Given F,(c,r) we wish to construct a Jeans distribution f,(m?, «) such that 


| | fasin 0, d0.db,=47 F,, 


"ful (resin @,)2, €) sin 0, d8, = F,(c, r). (9) 


Method.—Since p is known we may solve Poisson’s equation for V,,,, and 
add V ,,, (if any) to obtain V(r) as a known function. Now c?/2—V(r)=e and 
we write cr? =y, x/2r?—V(r)=e solving for r we obtain r=r(y,«), evidently 

C= on Se > 
r(x, €)}? 
note c being a magnitude >o. Write 


Fiic,r} =F, 1X —, r{x,<)}=Fly,e) 
© * (r(x) f 


Write ysin? @,= 7 then (g) becomes 
[’ f(T, «) dT = Fly,«) 
J 02xV1-T/x 
aCe 
Jo Vx-T 


* Or construct (following Kurth, Zs. f. Astrophys., 26, 168, 1949) the distribution function, f(«) 
dependent on energy only, which corresponds to the given density. 


=2V xF(x,€). 
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This is Abel’s integral equation and has solution 


f(Tye)= + [GN AE 
mJ o VT-x 
where the differentiation is performed at constant «. Thus we have found the 
Jeans distribution as required. We may remark however that most such 
distributions will be imaginary since we have no proof that f, is nowhere negative. 
We will have to throw away all these imaginary solutions and leave only the real 
ones before passing to the next step in obtaining the general distribution. 





Clare College, 
Cambridge : 
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PHOTOMETRY IN THE MAGELLANIC CLOUDS, II. SOME 
COLOUR-MAGNITUDE ARRAYS 


R. v. d. R. Woolley 
(Communicated by the Astronomer Royal) 
(Received 1959 September 16) 


Summary 


Colour-magnitude arrays are given for two small areas, and for some 
clusters, in the Large Magellanic Cloud. They are the result of measure- 
ment of photographic plates secured with the Radcliffe 74-inch reflector by 
Sandage and with the Victoria telescope of the Royal Observatory at the Cape 
by Alexander, and are referred to photoelectric sequences determined by 


Eggen and Sandage. The arrays show remarkable numbers of very blue 
stars. 





1.'Standard sequences in two areas in the Large Magellanic Cloud were 
determined by Eggen and Sandage and were described in Paper I of this series*. 
A number of photographic plates were exposed by Sandage with the 74-inch 
reflector at the Radcliffe Observatory, Pretoria, and by Alexander with the 
Victoria telescope at the Royal Observatory at the Cape of Good Hope. The 
latter is a twin telescope consisting of a photographic refractor of aperture 24 
inches and a photovisual refractor of aperture 18 inches, both having a focal 
length of approximately 23 feet. ‘The photoelectric sequences are used to cali- 
brate the magnitudes measured on these plates with the Sartorius Iris Diaphragm 
photometer of the Royal Greenwich Observatory. Generally speaking, the 
Pretoria plates are used to construct colour-magnitude arrays and the Cape plates 
to investigate variable stars, but some overlap of purpose has been necessary. 
The present note discusses four colour-magnitude arrays, as follows. 


(i) A field of sturs, including most of the photoelectric sequence LMC Field 
No. 1 of Paper I and measuring about 16’ x 14’ centred on R.A. 5" 00™, Dec. 
— 65° 34’ (1950), was examined. ‘The magnitudes of 327 stars were measured 
on a pair of Pretoria plates, the telescope having been diaphragmed to 43 inches. 
The exposures were 15 minutes both on the visual plate (103aD + GGr1 filter) 
and the blue plate (10320 + GGr13 filter). ‘The magnitudes found in this way, 
added to those of 21 stars in the photoelectric sequence, are plotted in Fig. 1, 
which shows B—V plotted against V. 


(ii) A similar field, of area 12’ x 11’ centred approximately on R.A. 5" 32™, 
Dec. —67° 08’ (1950), was chosen for examination. It includes some of the 
stars in LMC Field II of Paper I. Magnitudes were formed from measures on 
three Cape visual plates (103aD+QOmag 301) and three Cape blue plates 
(10320 +GG13). The results are plotted in Fig. 2. 


* Eggen and Sandage, M.N., 120, 79, 1960. 
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Fic. 1.—Colour-magnitude array, Field 1. 


(iii) The cluster NGC 2004 (R.A. 5" 29™, Dec. —67° 19’ (1950)) was 
examined. This cluster is described in NGC as “ globular, bright, pretty large, 
pretty rich, compressed, and having stars of twelfth magnitude ”. It is illustrated 
in Plate 6(a). This cluster is on the Cape plates of Field II and on four 
available Radcliffe plates. Sixteen stars marked in the illustration were chosen 
to form a sub-sequence and their magnitudes adopted by comparison with some 
stars of the Field II p.e. sequence occurring on each of two B and two V Radcliffe 
plates (the telescope being stopped to 43 inches). A total of 91 stars were mea- 
sured on three pairs of Cape plates and the two pairs of Radcliffe plates, relative 
to the sub-sequence as standards. The adopted magnitudes of the sub-sequence 
stars are shown in Table I. Of the stars measured, one proved to be an eclipsing 
variable. The resulting colour-magnitude array is shown in Fig. 3. 


(iv) Lastly, the neighbouring globular clusters NGC 1818 and NGC 1810 
were examined. NGC 1818 is described as “ globular, very bright, pretty large, 
rich, very much compressed and partially resolved”. It is situated at R.A. 
5" 04™, Dec. —66° 30’ (1950), and is illustrated in Plate 6(b). The nearby 
cluster NGC 1810 is “ considerably faint, small, round, and little brighter in the 
middle ’’. 
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Fic. 2.—-Colour-magnitude array, Field 11. 


Tape I 


Sequence for NGC 2004: provisional values 
B V 


m m 
11°70 12°06 


re g 
4 


15‘14 15°48 
13°38 13°66 
14°67 14°92 
11°75 12'10 
13°50 13°68 
15°10 15°34 
14°24 14°34 
14°88 15°17 
14°79 15°06 
14°56 14°86 
15°17 15°39 
15°36 15°68 
15°72 1§°70 
15°56 15°62 
15°55 15°74 
16°58 ~ 


QO 


vosg-r— rR Ao 
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. 3.—Colour-magnitude array, NGC 2004. 
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Fic. 4.—~Colour-magmtude arrays, NGC 1818 and NGC 1810. 
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These clusters appear on Radcliffe plates which have no p.e. sequence stars 
on them, and on Cape plates of Field I, which of course have the p.e. sequence, 
which is distant about 3} inches on the plate—the clusters being about 1} inches 
from the plate centre, and the sequence, in the mean, about 2 inches. (The 
exposed area on the plates is just under 6 inches square.) ‘The sub-sequence 
shown reproduced in Plate 6(6) and described in Table II was formed by taking 
the mean of about 10 transfers from the p.e. sequence on Cape plates. It is 
intended to check this sequence by setting up a fresh p.e. sequence near NGC 
1818 which may result in systematic corrections being applied to the results 
found for NGC 1818 and NGC 1810, but it is not thought that these corrections 
will exceed o™-10. The mean magnitudes found from four pairs of Radcliffe 
plates are shown in Fig. 4, those for NGC 1818 being shown as crosses and 
those for NGC 1810 being shown as circles. 


TasLe II 
Sequence for NGC 1818 and NGC 1810: provisional values 
Star B V B-V 
m m m 
13°35 13°28 + 0°07 
16°16 16°06 +o'10 
15°76 15°69 + 0°07 
15°50 15°36 +O°14 
14°92 14°96 — 0°04 
15°44 15°40 + 0°04 
16°04 15°86 +o'18 
16°32 16:05 +0°27 
16°57 16°52 + 0°05 


The magnitudes will be published im extenso on another occasion. In all 
cases, the results below visual magnitude 16 and photographic magnitude 16} 
must be regarded as rather uncertain—because of nearness to the plate limits and 
paucity of stars in the photoelectric sequences. However, Fig. 1 depends on two 
Radcliffe plates measured close to the p.e. sequence, and is probably reliable to 
another half magnitude. All the figures show a line indicating a (rather con- 
servative) estimate of the plate limit. 

2. All four fields show a remarkable sequence of blue giant stars and a few 
red giant stars, and a few stars of intermediate colour which are probably not for 
the most part members of the LMC at all. ‘The numbers of these stars of inter- 
mediate colour agree generally with the numbers of galactic stars per unit area 
given in Seares’ tables. 

The arrays of the two fields, and that of NGC 2004, are practically identical. 
The arrays of NGC 1818 and of 1810 are indistinguishable, and are displaced 
about o™-25 in colour and o™-60 in magnitude from those of the other fields, 
which suggests that there is slightly more absorption in front of these clusters 
than there is in front of NGC 2004 and the two fields. From the blueness of the 
stars in these three cases, the amount of absorption between us and these parts 
of the LMC must be supposed to be very small. 

3. The blue giant sequences, and the red giants at apparent magnitudes 13 to 
14 and colours about B—V=1™-8 resemble Sandage’s schematic array for h 


and x Persei. In terms of contemporary evolutionary theory, the areas and 
objects concerned are extremely young. 
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(a) NGC 2004: provisional magnitude sequence 





(b) NGC 1818: provisional magnitude sequence 
R. v. d. R. Woolley, Ph y in the Magellanic Clouds: 
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4. Some interest attaches to the question whether NGC 1818 and NGC 2004 
are globular. Since most, if not all, galactic globular clusters show red-giant 
luminosity arrays, as does for example M67, and are classed as “ very old” or 
“extreme population I” objects, some may choose to classify clusters by their 
colour-luminosity arrays. ‘To the present writer, however, it seems necessary to 
classify clusters (in principle) by referring to their internal motions, and to regard 
a cluster as globular if the great majority of the stars in it have not got the 
velocity of escape from the gravitational field of the cluster—and as open if the 
cluster is a moving association of stars, most of which have the velocity of escape 
from the others. If one accepts this dynamical concept of the difference between 
globular clusters and open clusters, and accepts the evolutionary view that clusters 
like M67 are red because they are old, one is led to the view that all globular 
clusters must have been blue once. Perhaps NGC 1818 and NGC 2004 are 
simply new globular clusters. ‘They have roughly the same diameter as NGC 
1783, which is a “ typical’ globular cluster without blue giant stars—namely, 
about one minute of arc, or about 18 parsecs—and the total brightness of NGC 
1818 appears (by inspection of the B plates) greater than that of NGC 1783. 
The latter may, however, contain the greater number of faint stars. 

It would be of great interest to determine the distribution of brightness of 
NGC 1818 with radius, to see whether it could be reconciled with the hypothesis 
that the stars in NGC 1818 are distributed subject to their mutual gravitation. 
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Summary 


Many investigations in radio astronomy are limited by the resolving 
power which can be achieved by conventional methods of aerial construc- 
tion. 

A new method of obtaining increased resolving power has been developed, 
which has been applied to the construction of both “ pencil-beam ”’ systems 
and interferometers. In this method two aerials are arranged so that their 
relative position may be altered to occupy successively the whole area of a 
much larger equivalent aerial. By combining mathematically the information 
derived from these different positions, it is possible to obtain a resolving 
power equal to that of the large equivalent aerial. Since the combination 
of the individual records may be done with different phase relationships, it is 
possible, without extra observations, to “‘ scan’”’ the synthesized aerial over 
an appreciable solid angle; because of this the total observing time of a 
synthesized instrument is of the same order as that of a conventional 
instrument. 

An interferometric system designed for the study of radio stars has been 
built which has an equivalent area for resolution of 8-105 sq. ft as well as a 
‘“‘pencil-beam’”’ system with an equivalent area of 3-:10%sq. ft. The sen- 
sitivity of both systems corresponds to a “‘ collecting area’’ of about 2-108 sq. ft. 





1. Introduction.—The operation of a radio telescope is limited by two 


dominating factors. One is the sensitivity, which specifies the smallest radio 


flux which may be detected, and the other is the angular resolving power. The 
relative importance of these factors is strongly dependent upon the type of 
investigation undertaken and for this reason specialized instruments are often 
developed for particular applications. In studies of the solar radio emission, 
for example, where high resolving power is needed at metre wave-lengths the 
primary requirement is that of angular resolution; the use of conventional 
aerials such as paraboloids is clearly impracticable. Even if structures large 
enough to give the required resolution were mechanically feasible their 
energy-collecting area would be far in excess of that needed to detect the relatively 
intense solar radio emission. Instruments have therefore been developed in 
which high angular resolution has been obtained by dispersing small aerials 
over a wide area of ground. Such systems include the diffraction grating type 
of array (Christiansen 1953; Blum, Boischot and Ginat 1958) and the 
interferometer of variable spacing (Stanier 1950, Machin 1951, O’Brien 1953); 
these two methods achieve the same result but in the former case all the aeriais 
are present simultaneously, while in the latter case two aerials only are moved, 
successively, to the different positions of the grating. 
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Similar techniques have been developed for the study of the distribution of 
the galactic emission. Conventional systems have been successfully employed 
for the short wave-lengths, but in many problems, particularly those involving 
the distribution of radio ‘‘ brightness ’’ at high galactic latitudes, it is important 
to make observations at wave-lengths of several metres; by dispersing the 
available collecting area it is again possible to obtain much larger resolving 
powers with adequate sensitivity. Two systems have been used, one of which 
employs a pair of mutually perpendicular thin apertures (Mills and Little 1953); 
the other utilizes a single thin aperture in conjunction with a small movable 
aerial (Blythe 1957). An earlier version of the latter instrument was used by 
Scheuer and Ryle (1953) to obtain increased resolving power in one coordinate 
for studying the distribution of radio emission across the galactic plane. 

The observation of radio stars presents a more complex problem in which 
three main approaches have been made. (a) The use of conventional systems 
of moderate size at a wave-length sufficiently short to provide the required 
angular resolution (Piddington and Trent 1956; Seeger, Westerhout and Conway 
1957; Roman and Yaplee 1958). These observations have been characterized 
by very small signal/noise ratios, and in consequence long integration times are 
used and extensive surveys have not been possible. (+) The use of conventional 
aerials of larger collecting area at longer wave-lengths (Hanbury Brown and 
Hazard 1953; Ryle and Hewish 1955); the signal/noise ratio of such systems 
has usually been large and the observations have been limited by resolution. 
(c) The use of ‘‘unfilled’’ systems such as the Mills Cross (Mills and Little 
1953) at a long wave-length. The angular resolution of such a system may be 
very large but the collecting area is smaller. 

If any major increase in the depth of radio star surveys is to be made it is 
necessary to construct systems with a collecting area exceeding 50000 ft® 
capable of operating at shorter wave-lengths or to use, at long wave-lengths, 
‘‘unfilled’’ systems having larger collecting areas than the Mills Cross. 

Since the flux density of the sources falls with wave-length and since the 
difficulties of constructing an aerial of given physical size get very much greater 
if it is to operate at short wave-lengths, it is possible that even with the new 
receiver techniques now being developed (Adler 1958, Giordmaine et al. 1958) 
it may be difficult to attain a sufficient improvement by using conventional 
systems at short wave-lengths. The alternative approach, which utilizes a long 
wave-length and correspondingly large collecting areas, appears to offer greater 
possibilities. 

In this paper an account is given of a new principle in the design of large radio 
telescopes. The method, which has already been used successfully by Blythe 
(1957), makes use of two aerials arranged successively in different configurations 
to provide information equivalent to that obtained from the use of an aerial of 
much greater physical size. Apart from a considerable economy of structure 
this method avoids some of the difficulties associated with the physical achievement 
of a graded excitation of amplitude and phase which is required in the case of 
large extended arrays such as the Mills Cross. Besides allowing greater collecting 
areas to be realized, the shape of the reception pattern can be adjusted, by 
computation alone, to suit different types of observation. The method necessarily 
involves considerable computation but this does not present a serious problem 
with the large electronic computers now available. 
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2. The principle of aperture synthesis.—It is well known that a simple 
interferometer, consisting of a pair of non-directional receiving elements 
connected to a receiver, gives a response proportional to one Fourier component 
of the two-dimensional distribution of radio brightness across the sky. By 
taking measurements in which both the spacing and the orientation of the 
interferometer are varied it is possible to determine the two-dimensional transform 
of the brightness distribution within certain limits and hence, by Fourier 
inversion, to derive the brightness distribution as observed with a specified 
resolving power (O’Brien 1953). Aperture synthesis may be regarded as a 
logical extension of this method but it will here be discussed in a somewhat 
different, but essentially equivalent, way which may give a more direct 
understanding of the method. 

Consider, for example, the operation of a receiving aperture such as a large 
two-dimensional array, or paraboloidal reflector. The signal which is delivered 
to a receiver from such an aerial may be regarded as the vector addition of the 
currents induced in each portion of the aperture. In the case of a paraboloid 
this vector addition is achieved, automatically, at the focal point, while for an 
array the same result is achieved by connecting all the elements to the receiver 
in the same phase. When the direction of an incident wave does not coincide 
with the normal to the aperture plane the induced current in the elements of the 
aperture will suffer a progressive phase shift and vector addition then gives rise 
to the usual directional properties of the aperture. In the case of an array, it 
it possible to add the currents with their phases suitably adjusted, and in this 
way to vary the direction of the principal response of the aperture; this feature 
is of great importance in large aerial systems where physical tilting of the aperture 
plane may be impossible. 

Now if it were possible to measure, separately, the current induced at each 
portion of the aperture by moving a small aerial successively across it, vector 
addition of the currents would give, for a constant source, exactly the same result 
as that obtained by using the complete large aperture. This process is not possible 
using a single receiving element since measurement of the phase of the induced 
current raises difficulties. By using two small aerials arranged as an interferometer, 
however, the relative phases are readily obtained and the successive determination 
of the current in this way forms the basis of the synthesis method. 

(a) The synthesis of a single aperture.—Consider a single rectangular aperture 
such as that shown in Fig. 1. The resultant current induced in the aperture by 
an incident plane wave is > /,,e*" where ¢, is the phase of the wave at the 


position of the nth element. The power P delivered to a receiver is then given 
by 


Px> Ene > ie o> 4, - = } Ae i cos (Din —, ). (1) 


If all the elements are of the same size, the first term in expression (1) is 
simply N times the power derived from a single element; the terms of type 
I,J, cos (¢,,—¢,,) are just the output obtained when a phase-switching receiver 
is connected to the elements m and n through equal cables (Ryle 1952). A single 
measurement of the power induced in one element, suitably combined with the 
outputs from a phase-switching receiver connected to a pair of elements, arranged 
successively to cover all possible combinations of positions, then gives a result 
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exactly equivalent to using the whole large aperture. By performing the addition 
with the value of (¢,,—¢,,) adjusted to give a progressive phase shift across the 
aperture it is possible to vary the direction of the maximum response of the 
aperture in a manner analogous to tilting the aperture plane. It is thus possible, 
without the necessity for further observations, to scan the synthesized aperture 
through an angle which is only limited by the directional properties of the small 
element. To perform the addition with modified values of (¢,,—¢,,) it is also 
necessary to measure terms of the type /,,/,, sin (¢,,—¢,,), and these are obtained 


by connecting the pair of elements in phase quadrature to the phase-switching 
receiver. 








Fic. 1.—The synthesis of a rectangular aperture. 


The measurements might proceed in the following way: having determined 
the power induced in a single element, an interferometer is set up with one 
element A fixed at position 1 and the sine and cosine terms are obtained by placing 
the other element B successively in all the remaining positions. Element A is 


then moved to position 2, etc., and the same procedure carried out. Now if this 
were done it is immediately obvious that many arrangements of A and B, such 
as those shaded in Fig. 1 are identical. A procedure for obtaining the different 
terms with no repetition is indicated in Fig. 2 (a) where the element A is fixed 
and element B is arranged to cover, once only, all the positions contained in a 
rectangle of approximately twice the area. It is of course possible to move both 
elements, in which case all arrangements can be obtained without moving the 
elements outside a rectangle of side D. It may be noted that if a total of N 


observations was required initially, the number now needed with no repetition 
is (2N)*2. 








Fic. 2.—(a) One method of synthesizing a rectangular aperture without repeating any element 
configuration. 


(b) The weighting function for each position of the movable element. 
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To obtain results equivalent to using the complete aperture the different 
terms must now be added with a weighting factor appropriate to the number of 
repetitions of each configuration of the elements. The form of the weighting 
function appropriate to each position of the movable element is sketched in 
Fig. 2(6). Direct addition of the terms without such a tapered weighting 
function gives a result equivalent to the use of an ‘‘ optimum’ array as discussed 
by Arsac (1955) and Barber (1958). While such systems give the greatest ratio 
of information to noise they are sometimes undesirable in practical measurements 
owing to the large subsidiary maxima inherent in their reception patterns. 

When the elements used in the synthesis have a width much greater than one 
wave-length, addition of the different terms with a progressive phase-shift is 
equivalent to the use of an aperture having a discontinuous “‘ stepped ’’ phase 
distribution across it. As in the case of an echelon diffraction grating this gives 
rise to more than one principal maximum in the reception pattern. This difficulty 
can be avoided if it is arranged that the different locations of the movable element 
overlap slightly. 

It is not necessary for the elements A and B to have the same shape and 
there are sometimes practical reasons why it is preferable for one of them to 
extend over the complete width of the required aperture so that synthesis is only 
performed in one dimension. A system of this type, which is exactly equivalent 
to that shown in Fig. 2 (a) is shown in Fig. 3 (a). When this method is adopted 
it is necessary to arrange that the current induced in each element of the extended 
array has the value appropriate to the weighting function shown in Fig. 2 (6). 
The aperture must consequently have a tapered distribution of excitation along 
its length. 














Fic. 3.—(a) One-dimensional synthesis of a single aperture. 
(6) One-dimensional synthesis of an interferometer. 


(b) The synthesis of an interferometer.—F or many investigations it is preferable 
to use an interferometer rather than a single receiving aperture and the method 
of synthesis may be applied in exactly the same way. By an obvious extension 
of the preceding arguments it is readily shown that an interferometer comprising 
two large apertures of side D may be synthesized as shown in Fig. 3(6). In 
the case of a phase-switching interferometer, however, there is no term 
corresponding to > /,,? in the output of the receiver and thus it is not necessary 
in the synthesis to measure the power received by a single element. 

3- Observation time and signal-to-noise ratio.—Before describing practical 
applications of the synthesis method it is necessary to consider how it compares 
with conventional methods in regard to the observation time and signal to noise 
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ratio. Since a synthesis using small elements requires many different measure- 
ments it might be expected that the observing time is necessarily longer. The 
time required for the repeated observations is, however, almost exactly 
compensated by the increased scanning rate which may be employed. It is 
useful to define a quantity which gives a direct comparison of the performance 
of synthesized and conventional instruments having the same angular resolution. 
In principle, the observation of a given area of sky can, by suitably arranging the 
scanning rate and the receiver time constant, be made to occupy a certain interval 
of time. We can therefore define the ‘‘ efficiency ’’ of a synthesized instrument 
as the ratio of its absolute sensitivity (i.e. signal/noise ratio) as compared with 
that of a conventional instrument of the same resolving power when the same 
observing time is allowed in each case. The efficiencies of various types of 
synthesized instrument will now be discussed. 

(a) Scanning rate and observation time.—During a survey a radio telescope 
is generally scanned continuously on one coordinate, either by steering the 
instrument or by using the Earth’s rotation, and a given region of the sky is 
investigated strip by strip. For our purpose it is more convenient to imagine 
that the scanning is performed in discrete intervals in both coordinates, the 
instrument being directed towards a certain point in the sky for the required 
integration time and them moved instantaneously to the next point. The angular 
distance between adjacent points in the scanning lattice depends upon the resolving 
power of the instrument and Bracewell (1956) has shown that the greatest interval 
allowed when the aerial is a rectangular aperture of sides a and 4 is given by A/2a 
and 4/26 in the respective coordinates, where A is the wave-length. No additional 
information is gained by using a finer scanning lattice but information is lost if 
a coarser lattice is adopted. 

Suppose, now, that a given solid angle © of the sky is scanned using a square 
aperture of side D. The total observing time is given by 4£2D*z/A*? where 7 is 
the integration time of the receiver. If the sky is scanned using the equivalent 
synthesized instrument in which the small elements are square apertures of side d, 
the time required to scan the solid angle for each of the arrangements is 422d?7/A? 
where it is assumed that the receiver time constant is the same as before. The 
number of different arrangements is approximately 2D?/d? and the total observing 
time for the same receiver time constant is therefore 8Q.D*7/A*. This is seen to 
be only twice that required had the complete aperture been used and it is important 
to notice that it is independent of the size of the small elements. 

To calculate the observing time when synthesis is performed in only one 
dimension, consider the system shown in Fig. 4(a@). If one of the elements is 
a rectangle of sides d, and d,, and the other an extended aperture of width d, 
and length 2D, then the scanning intervals are A/2d, and A/2D*. The total number 
of arrangements is approximately D/d, and the observing time is therefore 
40QD*7/\?. The observing time is thus half that needed when both elements 
were small and the same as if a complete aperture had been used. It is also of 
interest to consider the system shown in Fig. 4(6) which is similar to an aerial 
constructed by Mills except that one arm of the cross is omitted. Although no 


* The extended aperture of length 2D, when used with the appropriate tapered excitation is 
equivalent, when scanning, to a uniform aperture of length D. More exactly, the scanning interval 
is A\/2D+d, but we assume D>d,. 


16 





226 M. Ryle and A. Hewish Vol. 120 


synthesis is required such a system yields the same (non-repetitive) information 
as do the synthesized apertures under discussion and it will be shown that it has 
approximately the same efficiency. We note here that the scanning interval 
is A/2D, as for the complete aperture, so that the observation time is again 
40D? 7/22. 

It should be noted that this analysis assumes that the bandwidth of the receiver 
is not too great and that the area of sky under investigation is large enough to 
contain more than one coarse scanning interval A/2d. ‘The effect of a finite 
bandwidth is to narrow the reception pattern of the individual elements at large 
spacings (Ryle and Vonberg 1948), and this may necessitate the use of a finer 
scanning lattice. In the case of a single aperture it may be shown that the effect 
is negligible provided that Af/f is appreciably less than d/D; for an interferometer 
the ratio Af/f must be less than the ratio of d to the maximum element spacing. 
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(b) 


Fic. 4.—(a) The synthesis of a rectangular aperture. (b) A similar system which may be 
compared to the Mills Cross. It can be regarded as a synthesis in which all configurations are 
present simultaneously. 


(b) The efficiency of synthesized systems.—To derive the efficiency of a 
synthesized aperture we note that the signal obtained with a complete aperture 
of side D is directly proportional to its area D*. If, for an integration time 7, 
the root mean square receiver noise is m then the signal/noise ratio is proportional 
to D?/n. 

When the aperture is synthesized using small elements of area d? the signal 
obtained by adding the 2D*/d? different observations with appropriate weights 
is proportional to D*. Since the receiver noise is unrelated for the different 
observations the resultant noise is 2"*nD/3d where weighted addition has 
been performed. The signal/noise ratio is thus given by 3Dd/2!2n. Since 
however, the observing time for a fixed receiver time constant was twice that 
required when using a conventional aperture, the efficiency is given by 3d/2D. 
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In the case where a one-dimensional synthesis is used, as in Fig. 4 (a), the 
signal obtained for each arrangement of A and B is proportional to 2(Dd?d,)"* 
since the equivalent area is given by the geometric mean of the two apertures 
(Ryle 1952). The number of arrangements is D/d, and the signal/noise ratio 
is D(3d,d,)"?/n. In this case the observing time was the same as with a conventional 
aperture and so the efficiency is (3d,d,)"*/D. For the system of Fig. 4(6), in 
which no synthesis is employed, the signal is D(2d,d,)'* and so the signal/noise 
ratio is D(2d,d,)'*/n and hence the efficiency is (2d,d,)'/*D. All the systems shown 
in Figs. 3 and 4 thus have a comparable efficiency and the maximum possible 
efficiency in a synthesized aperture is (d/D)'* which is attained when d,=D. 

While the efficiency of a synthesized aperture is necessarily less than that 
of a conventional aperture of the same resolving power this should not be regarded 
as an immediate disadvantage. As mentioned earlier, the efficiency of a 
conventional aperture is often far higher than is required for a particular 
investigation so that the same physical aperture, arranged as‘a synthesized 
system of adequate efficiency and increased angular resolution, would have 
provided a much more powerful instrument. 

4. The practical application of the method 

(a) General considerations.—Since systems derived by one- or two-dimensional 
synthesis have approximately the same performance there might seem to be 
considerable advantages with a two-dimensional system arising from the economy 
of physical structure. There are, however, a number of practical disadvantages 
with two-dimensional synthesis and these are discussed below. 

(i) Data handling and computation: for a one-dimensional synthesis 
involving N arrangements of the elements, the corresponding two-dimensional 
system would require N? arrangements. The computational and data-handling 
problems are thus enhanced. 

(ii) Scanning: the angular scanning rate of a two-dimensional system must 
be considerably higher than for a one-dimensional system of the same resolving 
power for a survey to occupy a given observing time (cf. Section 3). It is a great 
practical convenience to employ the Earth’s rotation for scanning in right 
ascension, but this fixed scanning rate may be too slow for the desired observational 
programme. Facility for pointing the elements in both coordinates is then 
necessary which raises practical difficulties which may be avoided in a 
one-dimensional system. 

(iii) Receiver linearity: in a two-dimensional synthesis each intense source 
will be received over a solid angle A?/d?, and it may be difficult to ensure a 
sufficiently linear overall response to allow observation of weak sources within 
this area. In the case of one-dimensional synthesis the area of sky affected in 
this way is reduced by the receptivity of the long aerial to A?/Dd. 

(iv) Bandwidth: a limitation particularly relevant to the synthesis of an 
interferometer is set by the condition d/D>5f/F (cf. Section 3). Since the 
overall spacing of the synthesized apertures will be several times larger than their 
individual aperture D, the maximum bandwidth permitted in a two-dimensional 
synthesis is necessarily smaller than for the corresponding one-dimensional 
system in which the long aerial runs East-West. 

(b) The design of systems utilizing the Earth’s rotatton.—For the reasons 
outlined above systems adopting one-dimensional synthesis have a number 
of practical advantages. In current applications long thin apertures running 


16* 
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East-West, with provision for rotation about the long axis, may be constructed 
relatively economically and so the remaining discussion will be restricted to the 
design of such systems; it is possible, however, that other problems might reach 
a simpler solution if two-dimensional methods were employed. 

At first sight the restriction of the design to what is essentially a transit 
instrument, in which the scanning rate is fixed by the Earth’s rotation, might 
be thought to lead either to inadequate sensitivity due to a too short integration 
time, or to an unnecessarily long observing time. While it is true that the 
observing time required to survey a given area of sky is determined entirely by 
the angular resolution of the system and the rotation speed of the Earth, it is, 
however, possible to adjust the efficiency so that this scanning rate is always the 
optimum for a particular application. It is assumed that the area of sky involved 
includes 24" in right ascension and an angle of at least A/d in declination. 
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Fic. 5.—Part of a map obtained with the 7-9 m “ pencil-beam’’ system, 
The contours are in units of 1000 °K. 





In Section 3 the efficiency of a one-dimensional system such as that in Fig. 3 
was shown to be proportional to the quantity (d,d,)'?/D. Once the efficiency 
has been determined, the relative values of d, and d, may be decided by 
considerations of structural economy provided that their product is constant. 

5. The design of two practical systems 

(a) A pencil beam aerial for a wave-length of 7-9m.—This instrument was 
designed to provide a high resolution survey of the general galactic background 
radiation at a long wave-length so that, in conjunction with other observations 
at shorter wave-lengths, reliable information would be available on the spectral 
distribution of the emission from many galactic features. 

The original instrument built by Blythe (1957) employed one-dimensional 
synthesis as in Fig. 4(a) with an East-West array of length D= 1200 ft and a 
movable element of length d,=one wave-length. The effective width d, was 
less than the wave-length so that a single series of observations provided a 
complete map of the accessible sky. However, the presence of ground reflections 
led to an uncertainty in the receptivity in declination and gave an uncertainty 
in the overall scaling of brightness at southerly declinations. 
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The long fixed element of the radio star interferometer. 
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In the new instrument, which was designed to give an angular resolution 
of o°-8 x o°-8, a fixed aerial of length 3300 ft and width 40 ft was used so that 
a well-defined receptivity in declination was achieved at declinations above 
—10°. This increase in width and the use of a movable element of length 100 ft 
gave an efficiency of 7 per cent. 

A preliminary survey of limited resolving power (0°-8 x 1°-6), in which 
spacings of the movable aerial up to 30A were used, has been carried out and 
part of the map is shown in Fig. 5. 

















20" 40" 


Fic. 6.—A _ section of the first map obtained with the radio star interferometer. The most 
intense source at 204 37™ has a flux density of 11 x 10728 w.m.~2 (c/s)~1. 


(b) A radio star interferometer on a wave-length of 1-7 m.—A system capable 
of carrying out a deeper survey of radio sources than those made hitherto 
Mills and Slee 1957, Edge et al. 1959) was required, and considerable 
increases of both resolution and sensitivity were necessary. In addition the 
design was influenced by the importance of obtaining the positions of sources 
with greater accuracy, particularly in declination, and by the necessity for 
obtaining good measurements of the angular diameter of the most intense sources. 

An interferometric system having an approximately symmetrical primary 
reception pattern was adopted, in which the resolution in right ascension was 
provided by a 1450 ft parabolic trough running East-West (as in Plate 7). 
The required efficiency (about 25 per cent) was decided from an extrapolation 
of the number versus flux density counts of previous surveys; the ratio of 
the width d, of both aerials, and the length d, of the movable aerial were governed 
(cf. Section 4) by consideration of the relative costs of the parabolic elements 
and of the railway track supporting the movable aerial (Ryle 1960). 

A section of the first map obtained with this instrument is shown in Fig. 6. 
Some observations to establish accurate positions of 64 of the most intense radio 
stars have already been described (Elsmore, Ryle and Leslie 1960). 


Mullard Radio Astronomy Observatory, 
Cavendish Laboratory, 
Cambridge: 
1959 August. 
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THE ABSORPTION COEFFICIENT OF A PLASMA 
AT RADIO FREQUENCIES 


P. A. G. Scheuer 
(Received 1959 July 8) 


Summary 


The theory of radiation due to free-free transitions by electrons in a 
plasma has been the subject of much work. Elwert has given an exact theory 
for binary electron-ion encounters but there has been doubt about the 
applicability of the formulae when an electron moves in the field of many ions 
at comparable distances; it is the aim of this paper to remove that doubt. 

The power radiated by an electron at frequency v is proportional to the 
mean square amplitude of the Fourier component of frequency v of its 
acceleration; as Cohen, Spitzer and Routly have pointed out in a different 
context, one obtains the correct value of this mean square amplitude by 
treating all encounters as independent binary encounters, provided that the 
collision parameters concerned are all less than the Debye length. In the 
radiation problem, one is only interested in frequencies greater than the 
plasma frequency, and the binary collisions contributing to such radiation 
all have collision parameters less than the Debye length. Consequently 
Elwert’s treatment always gives the correct result, except at frequencies just 
above the plasma frequency. An explicit formula for the absorption coefficient 
of a plasma is worked out for a Maxwellian distribution of electron speeds. 

In the last section, the earlier considerations are applied to find the 
conditions necessary for radio-frequency radiation from encounters between 
electrons and concentrations of space charge. 





1. Introduction.—The emission and absorption of radio waves by ionized 
hydrogen is of considerable importance in astronomy, both in connection with 
the radio emission from the solar corona and in connection with the emission 
(especially at centimetre and decimetre wave-lengths) and absorption (at 
wave-lengths of 10m and more) of radiation by interstellar H11 clouds. With 
the recent development of low-noise amplifiers, it may be expected that greater 
precision will be achieved in high-frequency observations than has hitherto 
been possible, and particularly for this reason it seems desirable to have an 
accurate theory of the emission and absorption of radio waves in a plasma. 

Much has been written on the subject, but I shall not attempt to give the 
history of the problem, and refer only to the work of Elwert (1) and of Smerd 
and Westfold (2). Elwert’s work deserves particular attention, because he 
derives an exact quantum mechanical formula for the radiation from collisions 
between electrons of initial velocity v and positive ions of charge Ze. He then 
shows that the exact formula contains Gaunt’s quantum mechanical formula (3) 
as one limiting case, and a classical formula as another limiting case; he also 
shows that the classical formula is an excellent approximation to the exact formula 
in radio-astronomical applications. Elwert’s result has, however, been open to 
one objection. Occasionally an electron moving through a plasma passes so 
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close to an ion that its motion is controlled predominantly by the field of that 
ion (and much of the radiation is emitted in these close encounters), but 
generally the electron moves in a field to which several ions make comparable 
contributions, so that a theory which considers only binary collisions is not 
applicable. Smerd and Westfold took the view that, since the fields of different 
ions will, in general, tend to cancel out, one must consider only encounters with 
collision parameters 5 less than the mean interionic distance (b6,=N;~"*). 
Collision parameters do not occur explicitly in Elwert’s quantum mechanical 
calculation, and Smerd and Westfold adopt a different approach, which is purely 
classical. They also introduce the effect of the refractive index of the plasma for 
the first time. Smerd and Westfold’s procedure is not convenient for the present 
discussion, so it is proposed to show in outline how the classical formula for the 
absorption coefficient may be derived from first principles, taking ‘‘ multiple 
collisions ’’ into account, and at the same time carrying out the calculation to 
a better approximation. 

2. Definitions.—Various quantities which will occur in the calculation are 
defined below, in particular, various characteristic lengths (b, to b,) of the plasma. 
v denotes the speed of an electron, —e its charge, m its mass. N, and N, are the 
numbers of ions and electrons, respectively, in unit volume. Ze is the charge 
of an ion; & is Boltzmann’s constant. T is the kinetic temperature of the 
electrons. 4 is the collision parameter, which is defined for a binary encounter 


as the distance of the ion from either asymptote of the electron’s (hyperbolic) 
trajectory. 


b, is a value of 6 chosen to be >b, and < various upper limits (see below). 


b, = N,-"8, the mean interionic distance. 


b,= ; the distance travelled by a typical electron in one period of the 


radiation (frequency v) under consideration, is 27},. 

mae 3kT U2 a (v?)'2 
47 Z*e* 279 
expression is the plasma frequency. 


b, is given by at 3RT= }m(v?). 
bb, 2 


is the Debye shielding length; v, in the second 


At distances <5, from an ion, an electron’s potential energy is comparable 
with its kinetic energy. An encounter with collision parameter <), scatters the 
electron through a large angle. 

3. Description of the physical processes—Suppose now that an electron is 
moving through a plasma; provided that the kinetic temperature of the ions is 
not much greater than that of the electrons, the electrons move much faster than 
the ions, and (to the accuracy to which we are working) the ions may be regarded 
as stationary. 

Every time an electron passes close to an ion, it is deflected; neglecting 
radiation damping, it performs a hyperbolic orbit (Fig. 1(@)). During this 
encounter, it suffers accelerations whose components vary with time as shown 
qualitatively in Fig. 1 (b) and (c). The frequency spectra of these acceleration 
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components are also illustrated; the ‘‘transverse’’ acceleration, Fig. 1 (6) and 
1 (d), has a white spectrum up to frequencies of the order of v/b but cuts off at 
higher frequencies. The ‘‘longitudinal’’ acceleration is only important at 
frequencies of the same order of magnitude as v/b. 











ae ee 


oa _roaey v ; — v 


Fic. 1. 





We shall want to calculate the radiation from the encounter. The electron, 
together with part of the charge on the stationary ion, constitutes an electric 
dipole of varying moment; this dipole may be regarded as a superposition of 
sinusoidally varying dipoles. If the moment is p=), cos 27vt, then p is simply 
the electronic charge times the appropriate Fourier component of the acceleration, 
and the power radiated in a medium of refractive index , is 


(p= — = (mean squared amplitude of p). (1) 


3H rr 

Substitution of the aj accelerations (Fig. 1 (b), (c), (d), (e) in (1)), 

followed by integration over all collision parameters, should, if correctly 

performed, yield Elwert’s classical approximation for a binary counter. However, 

the need to consider multiple collisions makes it more convenient to split the 
problem into two parts. 


(i) ‘‘ Distant encounters’’, including multiple encounters 

In the sort of plasma we shall consider, the kinetic energy of a typical electron 
is very large compared to its mean negative potential energy in the field of an 
ion (i.e. 6,>5,), so that the electrons travel approximately in straight lines, 
with occasional sudden changes in direction due to close encounters with ions*. 

So long as the electron moves almost uniformly in a straight line, the 
accelerations due to the electrostatic fields of the ionst of the plasma may be 
considered as small perturbations on a linear path. The acceleration experienced 
by the electron is then, to a good approximation, 


e 
— E(t), 
= E(t) 
where E(t) is the field which would be experienced at time ¢ by an electron moving 


through the plasma with uniform velocity. 


* The condition 6, > 6, requires N;~!*> Fo that is T>10~*ZN,' thus for hydrogen at 10° K 


(or hotter), 6,/b,2107N,~! and the condition is valid by a large margin up to N, =10"* cm~ 
+ Electron-electron collisions do not normally contribute to the radiation; see Section 7. 
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Using this approximation, the ‘‘ multiple collisions’’ can be dealt with. 
We shall Fourier-analyse the field E(t) due to a random distribution of ions, 
ignoring those ions which lie within a distance 5, (6, >6,>6,) of the electron’s 
path; the value of p corresponding to a Fourier component of E(t) is obtained 
by multiplying by e.(e/m), and can then be substituted in (1) to obtain the 
power radiated. 

When an electron approaches an ion within a distance of the order of 4,, 
it is deflected through an appreciable angle, and the power contributed by the 
pulses of radiation emitted during close collisions must therefore be calculated 
separately and added. 


(ii) ‘* Close encounters”’ 

The pulse of radiation emitted during a close encounter lasts for a time of 
the order of b/v or less. If b/v<1/v, v lies on the white part of the spectrum, 
Fig. 1 (d), and one can approximate to the actual pulse by an infinitely short pulse 
whose magnitude is determined by the total momentum change suffered by the 
electron during the encounter. Also, with b/v<1/v the contribution of the 
longitudinal accelerations is negligible (Fig. 1 (e)). 

Now the condition 6/v<1/v is equivalent to b<b,. If b,<b,, then almost 
all the encounters capable of emitting radiation of frequency v are binary 
encounters, so that a further cut-off at a collision parameter 5, need not be 
considered, and Elwert’s formulae are clearly valid. If 5, is greater than, or 
comparable with, b,, then b <4, implies <6, ; in all the cases we need to consider, 
b<b,<b,, and hence b<6,. ‘Thus we can calculate the contribution of the close 
encounters (b <b,) using the assumptions that 


(a) the encounters are binary, since b<b,; 


(6) the pulses are infinitely short, since b<4,. 


4. Calculations.—-We now calculate the radiation from the two classes of 
encounters separately, dividing the collision parameters into ranges above and 
below 5). We shall then integrate over a Maxwellian distribution of speeds. 

(i) Distant encounters—b > by 

Consider an electron moving along the x axis, and an ion at x=o0, y=6, z=o. 
The electrostatic potential V due to the ion at a point on the x axis is (4) 


V = Ze(b* + x*)-!? = Ze | : Khy(£4) cos &x dé. 


0 


In the second expression, V’ has been Fourier-analysed into cosine components 
along the x axis with wave-numbers £. Kh is a Bessel function of imaginary 
argument, of the second kind, in the form used by Heaviside and by Jeffreys. 
The field components are therefore 
oV 5 


Se oe wit | " & Khy(€b) sin Ex dé, 


Ox 
ss oe = Ze | “€ Kh,(£b) cos éx dé, 
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Suppose we consider only the Fourier components of the field with wave- 
numbers & to +6. Then 


E, (€ to £ +8) = Ze } re Khg(€b) sin éx dé 
7 € 


sin 45éx 


= Zeé Kh,(£5) sin Ex dg : bdEx 


and similarly 


E, (£ to € + 8¢) = Zeé Kh, (£b) cos £x 8. ee 


A simple transformation of coordinates shows that for an ion at (X, bcos ¢, bsin¢) 
instead of (0,b,0), we should obtain 


sin 4dé(x — X) , , 
E, (€ to € + 3€) = ZeE Khg( Eb) “= = 8€ singe —X), 


on ele cos ¢ 5€ cosé(x— X), 
sin | bdE(x — X) X) 
~ $8E(x =X) 
The field actually experienced by the, electron is the sum of fields FE, due to 
many ions; these lie at various distances X; along the x axis, at various distances 
b, from it, and various angles ¢; around it. We shall require the mean square 

of each component of the resultant field, for example 


{> E,; (€ to £+d£)P = YE AE to €+5£)+ 23 E,, (€ to €+d€)E,, (€ to & +98), 
If the positions of the ions are uncorrelated, the field components due to the ith 
and jth ions are equally likely to have the same or opposite sign, and therefore 
the mean value of } E,,E,, is zero. The mean number of ions in the typical 
volume element is N,.db.bdd.dX, and hence 


E,, (€ to €+0£))?= SE,P 


E, (€ to + 6€) = Zeé Kh, (£b) 


E, (€ to £ + 5£) = Zeé Kh, (£6) — sin 6 5€ cos (x — X). 











=N,|° Z%%€* Kh,2(¢b)b db 3 cos? pd 


be 


x “s sin® 46€(x — x) Os -_ 
| w UBE(x— xy cos*® £(« — X) dX (5£)* 


= ZN, in Kh,*(£b)(E6) d(€b).. § 755. (5£)?, since 5¢ <¢ 
bo€ 


= ZNSE [- 7 Kh,*(y) dy (2) 
bo € 


and similarly 





(SE, € to € + 0€)}* = 20°Z*e2N, 5£ - n Kh,2(n) dn. (3) 


The mean squared z component is equal to the mean squared y component. 
The component of the field with wave-number £ causes periodic acceleration 
of the electron with frequency v={v/27; the magnitude of the acceleration is 


< Eq (€ to € +66) 
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in the x direction, etc., and hence the mean power emitted according to (1) is 
_2_ F+5,3+5 - 2,{(£) SE Eyi+... | 
pa OP+BF+b2) = 1(5) (SEs ewer sey +...| 
2e* ownanr [" ; : 
se 722 f 2 2 

by (2) and (3) 

"s 4n*Z*eAN S(2mv/v) (” 


3Hem* 
Integrating by parts, 








7{ Kho*(n) + Kh,?(n)} dy. 


oe by 2aviv 


[7 Kb.%(n) dn = [Khy(n)(—7 Khg(n))]— { (~n Kha(n))( ~ Kbg(n)) dy 
that is, 


[ n(Kh,2(n) + Kh,*(n)} dy = [— Khy(n) Kh,()}. 


_ 
Since Kho(y) and Kh,(») both decrease exponentially as 7 oo we obtain: 
mean power radiated by an electron due to encounters with b> 6, 


F26N' 5 , , 
i ae (72 bn) Kh, C= by) Kh, (22 by) (4) 
3p mit c i t 


(ii) Close encounters—b < by, 


A binary encounter with collision parameter 5 deviates the path of an electron 
through an angle 26 such that 


sin 6 = {1 + (bmv*/Ze*)?\~)2, 
Hence the electron suffers a change in velocity of 2vsin@, and 
Ap= | Bae =e.2usin 6 = 2ev{1 + (bmv?/Ze*)*}-!2, 
We now have to Fourier-analyse a sharp pulse in p, the magnitude of the pulse 
being given by the above value of | pdt. Isolating the components in the 
frequency band v to v+éy, we have 


v+ by 
p (v to v+dv)=2Ap [ Cos 27vt dv 


= 2A p dv (sin rdvt/mdvt) cos 2n(v + 45v)t. 
This represents an oscillation of frequency v+45v with a slowly varying 
amplitude 2Apsin wévt/mt. By (1), the total energy radiated is 
Te ts sin 7b vt\? I ; 
—-, dt = — 2 8p. 
5a | ~« (23 mt ) 3 4OP) ” 


An electron moving with speed v has 27bdbvN, encounters per second with 
collision parameters 6 to b+db; hence the close encounters (b<5,) made by 
an electron with speed v contribute a mean power 


» b 25,2 
a a oe e 4e°t a 
. oa (Ap)*27b db uN, dv = [ ie Ty ome zea 2mbuN,db dr 


167Z7e8N mv*b,\? 
= ———‘In({ 1+ (— bv 
3pemie Ze 





in the frequency band » to v+ dv. 
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(iii) Integration over a Maxwellian distribution of electron speeds 
The next step is to integrate both (4) and (5) over a Maxwellian distribution 


of speeds. For this purpose, neither expression is convenient, and the assumption 
b,>b,>6, must be used again to simplify them. In (4) 


and for all but the slowest electrons (about (4,/b,)° of the total) this number 
is <1. Therefore one can approximate to the Bessel function in (4) by 


, 2 , 
Khy(n)=—<(Indyty);  Khy(y)== — (y=0'577...) 


TT 
obtaining : 
(LON S v 
Mean power = pee In —y)}. 
zu mev trvb, 
In (5), 
mu?, 20° 5 
Zéaé* oo b 
which is >1 for all but the slowest electrons (about (4,/b,)°? of them), so that 
one may replace (5) by : 


Zie6 N. 2 
Mean power = ie In (a b ) ‘ (5a) 
Adding (44) and (§ a) and integrating over a Maxwellian distribution of speeds, 
we obtain the power radiated by the electrons in unit volume of the plasma: 
2kT 


zum nZev 


Using the formula 


0 Uv 


F2B®N Sy (1 [ 3 3/2 mrt 
320Z ca | = {In me r} 4 oy ( m ) vte~ BT dv. 


— WN, 
vr 


| “Inne-"dn= —y, 

0 

the integration can be performed, giving 
32nZ*e®N; N Sv ‘in 8k°T3 - (6 
3u0m(anmk iin} mieim >" 


To obtain the corresponding absorption coefficient, A, we note that the power 
radiated from unit volume is K.87k7 *5v/c?, hence 


8nZ*eN;N, { 8k T? } 
In —Sy>. 








~ 3uc(ammkT P22 || 22m (7) 


The ‘‘ logarithmic factor ’”’ 


tin moe } 

nZey 4 
which occurs when (44) and (5 a) are added to obtain the mean power radiated 
by an electron is precisely the same as the factor occurring in Elwert’s formula 
for the effective cross-section ((1), end of p. 480, formula for g’), and (7) is the 


formula which is obtained on integrating Elwert’s results over a Maxwellian 
distribution. 


On comparing (7) with Smerd and Westfold’s formula 


- 8a Z*tNN, 4kT, \? 
K~ Jdark TPs (: +(44)) 8) 
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one finds that the essential difference is that in (7) 5, has been replaced by dy. 
That is to say, the upper limit to collision parameters does not occur at the mean 
interionic distance, but when the pulse of radiation becomes so long (> 6b,/v) 
that it contains no appreciable Fourier component of frequency v. Since (7) 
agrees precisely with the classical approximation to Elwert’s binary collision 
formula, it appears that even in the multiple collisions, with b> 5,, the interaction 
between the electron and each ion contributes to the radiation as if it were a 
separate binary collision. This last fact may well be regarded as obvious if one 
considers that the radiation is proportional to p?, and hence to the square of the 
field to which the electron is exposed, and that the mean square resultant of a 
number of vectors (the fields of the ions) is equal to the mean sum of the squares 
of the individual vectors (cf. Cohen, Spitzer and Routly (5) who make the same 
point in connection with electrical conductivity). The above calculation has 
nevertheless been set out in some detail, because the matter appears to have been 
in doubt for some time, and also because the formula (7) does not appear to have 
been published explicitly. 

5. The numerical values of K given by the formula (7) do not differ drastically 
from those given by Smerd and Westfold’s formula; for conditions of 
astrophysical interest the difference rarely exceeds 25 per cent. Perhaps the most 
interesting feature of (7) is that it predicts a slow decrease in the flux densities 
of thermal radio sources at high frequencies, rather than the perfectly white 
spectrum given by (8). For a nebula of protons and electrons at 10°K and 
over the frequency range 300 Mc/s to 3000 Mc/s, the change would correspond 
to a spectral index (6) of —o-1. 

Owing to the large overlap between the regions of validity of the approximations 
used for close and distant encounters, formula (7) is a fairly precise one, and 
should be adequate for foreseeable radio-astronomical applications. The largest 
error is the neglect of Debye shielding (see Section 6); this introduces a 
fractional error of the order of v»/vIn(6,/b,), which might amount to a few 
per cent near the turning point of a ray in the solar corona. 

6. The crucial assumption in the calculation is the assumption that the 
positions of the ions are uncorrelated. This assumption requires critical 
examination, for the shielding of ions (originally investigated in the case of 
strong electrolytes) and the existence of collective phenomena show that it is, 
for many purposes, quite untrue. 

The situation indicated by (6), that the mean power radiated is the same as 
if all electron-ion encounters were independent, may be pictured as follows: 
when an electron moves in the fields of several ions, these fields, in general, tend 
to cancel out, so that the radiation per encounter is reduced. There are, however, 
statistical fluctuations in charge density, and a region containing an excess of n ions 
will act like an ion of charge nZe on electrons passing near it. The power radiated 
is proportional to the square of the ionic charge, so that the cluster is m times as 
effective as the n ions individually. As shown in the previous section, the two 
effects cancel out exactly when the distribution of ions is random. In a real 
plasma, the probability that a particular kind of cluster of ions will occur is the 
product of its statistical weight and a factor exp ( — «/k7T) where € is the electrostatic 
energy required to form the cluster. The average number of ions in a cube 
of side r is r°N,, and this number is subject to fluctuations of about (r°N;)"?; 
fluctuations greater than 5(r°N;)'* are very rare for statistical reasons alone. 
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The fluctuations in a real plasma are appreciably smaller than they would be 
in a random plasma if and only if the electrostatic energy 


{Ze(r9Nj)!2}2 
r 


«= } (charge)* (capacity)~' = = Z*e*7? N; 


is greater than or comparable with k7, that is, if 


r>kT/ZeN,. 
But 


me ( 3kT :) er (v2)? 
4nN,Z7%e? 27) 

is the Debye shielding length; therefore we conclude that fluctuations on a 

smaller scale than the Debye length are unaffected by electrostatic forces; 

fluctuations on a larger scale than the Debye length are much smaller than they 

would be in a random distribution of ions. (Pines and Bohm (7) give a thorough 

discussion of fluctuations in a plasma.) 

Since radio noise can only be emitted if the refractive index p is real, the 
frequency v radiated is always greater than the plasma frequency vg, and therefore 
b, = (v?)"2/27v is always smaller than 6, = (v*)"2/27v,._ But a charge concentration 
is only relevant to our problem if its diameter is of the order of 5, or less; an 
electron would take longer than 1/v to pass by a larger charge concentration. 
The assumption of a random distribution of ions is therefore correct so far as 
the emission of radiation is concerned, with a possible error of the order of one 
in the logarithmic factor of (7) when v is only a little greater than the plasma 
frequency vp. 

Two minor assumptions which are implicit in the calculation may also be 
dealt with here. Firstly, it was assumed that the pulses of radiation, due to close 
encounters, occur at random intervals, and that the energies radiated by them 
may therefore be added. A collision with b<b, may be expected to occur once 
in (7b,?.N;)~1 = (6,/by)?(b,/7) of electron trajectory. The exact positions of two 
ions so far apart could not be affected by their electrostatic repulsion, since their 
mutual potential energy Z7e?(b,/by)~*(by/7)~! = Z4kT2(b,/b, )?(6,/6,) is much less 
than their thermal energy. ‘To affect the relative phases of the radiation on 
frequency v from the two pulses, ordering of the ions must have a precision of 
the order of 5,, and, since b,<;, the over-dispersion on a scale >b;, or even 
collective phenomena due to non-thermal disturbances, are irrelevant. Secondly, 
it has been assumed, in calculating the radiation from distant encounters, that the 
electron’s approximately linear trajectory continues indefinitely. In fact, the 
trajectory, and hence the wave-train radiated, is broken up into sections by 
close collisions. A rough estimate of this effect may be obtained by regarding 
these lengths of wave-train as uncorrelated in phase, as in the simple theory 
of the collision broadening of spectral lines; we then find that the spectrum of 
the radiation is smoothed over a frequency range of about 


0(b,/by)-*(by/7)* = 277? v(bo/b,)?(b2/b,), 
which is a very small fraction of v. The unsmoothed spectrum, given by (6) 
contains no sharp resonances and the effect of the smoothing is therefore 
negligible. 
7. Radiation due to large concentrations of charge.—The considerations of 
Section 6 lead naturally to a few comments on the possibility of radiation by 
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encounters between electrons and charge concentrations induced by non-thermal 
disturbances. ‘This possibility was considered by Hoyle (8). 

It was shown in Section 6 that, in a plasma in thermal equilibrium, statistical 
fluctuations in charge density are normal over regions less than 4, in diameter, 
but fluctuations on a scale greater than 6, are much less than in a random 
distribution of ions in space. Conversely, scales greater than the Debye length 
are the domain of collective phenomena; plasma oscillations and other collective 
phenomena on a scale smaller than the Debye length would be quickly damped 
out by diffusion. Now suppose that charge concentrations have been set up in 
a plasma by some external disturbance: these charge concentrations will be very 
much more effective in scattering electrons than single ions would be, owing to 
their greater charge. A charge concentration is at least b, in diameter, so that an 
electron takes a time > b,/v to pass by it, and might be expected to radiate a pulse 
of about this length. However, b,/v is of the order of (27v,)~! for a typical 
electron of the plasma, so that the ‘‘pulse’’ contains no frequencies which 
can be propagated. Hoyle abandoned his suggestion for this reason. Radiation 
could only be emitted by means of a ‘‘ space-charge concentration mechanism ”’ 
if some electrons much faster than the thermal electrons were present. 

It is not inconceivable that favourable conditions could occur for short periods 
in nature. In particular, Kahn (g) and Buneman (10) have shown theoretically 
that when two plasmas collide, the energy of mass motion is converted into 
chaotic particle motions via plasma oscillations of large amplitude and on a scale 
only a little larger than the Debye length. Complete charge separation on the 
scale of the Debye length is the best possible condition for radiation by the 
space-charge concentration mechanism. Buneman has even suggested that 
electrons with more than thermal] velocities may be generated in the same process. 
It is also known from observation that streams of fast particles occur in various solar 
disturbances. 

In estimating the radiation emitted (or the corresponding value of the 
absorption coefficient) one need only consider encounters with collision parameters 
greater than the radius 7 of a typical charge concentration, for an electron passing 
through the charge concentration does not experience the full Coulomb force 
corresponding to the total charge concentrated at the centre of the charge 
concentration. Also, the effective b, is always less than b,, and r>6s, so that 
the approximate method (i) of Section 3 can be used; putting 6,=r in (4a) 
one obtains the mean power per fast electron (v is now the speed of the fast 
electrons, Ze the total charge of a charge concentration and N, is the number 
of charge concentrations per unit volume (about }7° if they fill space)). The 
formula for the absorption coefficient, analogous to (7), is easily derived; it 
differs from (7) only in the interpretation of 7, Ze, and N,, and in the logarithmic 
factor, which is essentially 2 In (b,/6,) in (7) and becomes 2 In (5,/r) for the present 
case (using the speed of the fast electrons in the definition of 5,). 

The spectrum of radiation emitted by a space-charge concentration mechanism 
would be broad, and could not contain harmonically related frequency bands, 
such as are observed in some types of solar burst. Also, if a frequency n times 
higher than the plasma frequency is to be radiated, electrons with at least n times 
the speed of the thermal electrons must be present. 

Normally, electron-electron encounters (and ion-ion encounters) can be 
neglected, for the two electrons suffer equal and opposite accelerations, and hence 
the dipole radiation from the two electrons interferes destructively. The 
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quadrupole radiation from the pair of electrons is important only if their 
separation 6 is comparable with the wave-length of the radiation; but the 
radiation emitted in the encounter extends up to a frequency w/b, and hence 
down to a wave-length of (c/v)b, which is >6 unless v is comparable with c. 
For relativistic electrons, the electron-electron encounters might contribute 
appreciably to the radiation. 

The absence of radiation from electron-electron encounters has one 
entertaining consequence. Suppose a fast electron moves in a plasma in which 
there are concentrations of ions, but that these concentrations also have enough 
electrons associated with them to make them substantially neutral. According 
to the above argument, the presence of the electrons should not affect the amount 
of radiation emitted by the fast electron, but common sense insists that the 
electron cannot be scattered significantly by a blob of neutral plasma. ‘The correct 
physical picture in this case must be that the fast electron is virtually undeviated, 
but it repels all the electrons in the blob of high density, and these electrons 
therefore radiate coherently. In other words, the fast electron temporarily 
polarizes the entire blob. Thus real charge concentrations are not necessary 
for radiation by the ‘‘ space-charge concentration mechanism’’; irregularities 
of density, such as might be produced in shock waves or sound waves, would 
be just as effective, if they could be formed on a sufficiently small scale. 

However, in astrophysical examples, the mean free path (the smallest scale 
of sound waves) is very much larger than the Debye length, so that this version 
of the mechanism is unlikely to be common in nature. If a magnetic field 
permeates the plasma, the smallest possible scale of “‘ sound ” waves is the Larmor 
radius, but even then a fairly large field (such as that over a sunspot) is needed 
to make ‘‘ useful’’ irregularities possible; in a field determined by equipartition 
between magnetic and thermal energy, the Larmor radius would be c¢/vp, the 
‘**plasma wave-length’’, which is much too large. 
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Summary 


A high resolution survey of the Cygnus X region was made at 408 Mc/s 
using the 250 ft paraboloid at Jodrell Bank. The results were combined 
with the high resolution surveys at 1390 Mc/s of Drake and Westerhout to 
investigate the spectral characteristics of the individual radio sources in 
Cygnus X. The results support the previous conclusions that the emission 
is from optically thin Ht regions except for the strongest source in the 
region, the y Cygni source, which was found to have a non-thermal com- 
ponent. 





Introduction.—All surveys of the radio emission of the Cygnus X region 
have indicated that the origin of the radiation is in an extended H11 region. 
Piddington and Minnett (1) first showed that the source had a thermal spectrum. 
Davies (2) has confirmed this theory and also concluded that the source is located 
in the 2nd spiral arm of the galaxy. Recent high resolution surveys at 1390 Mc/s 
by Drake (3) and Westerhout (4) have resolved the region into a large number of 
sources. 

The 250 ft paraboloid at Jodrell Bank (5) was used to survey the Cygnus X 
region at a frequency of 408 Mc/s. As the beam size was the same as that used 
by Drake at 1390 Mc/s, a detailed investigation of the spectral characteristics 
of the individual sources was possible. 


1. Theory 


(a) Flux measurements at 408 Mc/s using the 250 ft paraboloid.—The termi- 
nology in this section is that given by Seeger, Westerhout and van de Hulst (6). 

The full beam pattern of the radio telescope was obtained using the Cassiopeia 
radio source, from which the effective solid angle 22’ of the beam was calculated 
to be 0°68 sq. degrees. ’ is given by the expression 


Q’= | £(0,4) dO 
full beam 

in which f(6, ¢) is the power received in the direction (8,4) normalized by putting 
f=1 atthe centre of the beam. The full beam is defined by a circle with diameter 
5 times the half-power diameter of the beam. 

The sensitivity was checked daily by measurement of the deflection D, 
produced by the 14N5A radio source. This source has a flux density S, of 
47 x 107% wm? (c/s)! (7). 
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The flux density S of any distributed region bounded by the limits («,, a) 
in right ascension and (8,, 5,) in declination is then given by 


Pay (es 
[of (8) daa 
Xd 4, 
Q2’Do 
where D (a, 5) is the deflection at position («, 4). 
From these flux measurements the brightness temperature 7, of a uniformly 
intense region is derived by means of the standard formula (8) 
<a 
2kQ2, 
where Q, is the solid angle subtended by the source. Thus the brightness tempera- 
ture is related to the deflection D by the expresssion 
4 SA? D 
TS 
2kQ’ Do 

(b) Spectral analysis.—In comparing surveys of regions at differént frequencies 
with telescopes of different size, the brightness temperature of an extended 
region depends only on the spectral characteristics of the emission. For point 
sources the brightness temperature depends not only on the spectral character- 
istics of the source but is proportional to the size and diffractive efficiency of 
the aerial. 

Seeger, Westerhout and van de Hulst (6) have shown that for a telescope 
of effective solid angle ©’ an increase of flux from a point source of 2kQ/A? at the 
centre of the beam is equivalent to a 1°K increase in 7), over the full beam. 

Two spectral classes of radio sources are commonly recognized, i.e., those 
with a thermal spectrum and those with a non-thermal spectrum. For thermal 
sources the radio brightness (flux density per steradian) has a spectral index 
which may range from + 2 for optically thick regions to o for optically thin regions. 
For non-thermal sources the spectral index usually lies between o and —1. 

A comparison of the spectral characteristics of the sources in Cygnus X is 
made in the following sections using the 20 cm survey of Drake (3) and the 
present 75 cm survey. The effective solid angle of the aerial beams was the same. 
Table I gives the ratios of brightness temperature to be expected in these two 
surveys for various classes of radio source. 


S=%,. 





Tasie I 
Spectral index 
Class of source Brightness 

flux/ster. ve 
: b 

density 

Extended optically 
thick thermal +2 


Extended optically 
thin thermal 


Extended non-thermal 


Point non-thermal 
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2. Results.—Fig. 1 shows the 408 Mc/s radio brightness contours of the 
Cygnus X region. Fig. 2 is a typical scan at constant declination + 43°45’. 
The sources are listed under their W esterhout numbers (4), together with positions 
and peak brightness temperatures, in Table III]. The criterion used in the 
recognition of a source is that it should be possible to draw a closed contour 
comparable with the beam width around the point. All the temperatures have 
been measured above a reference region (R.A. 20"58™, Dec. + 39°) which is 
5° from the galactic plane. 





Declination 











Right Ascension 


Fic. 1.—408 Mc/s radio contours of the Cygnus X region, indicating the brightness temperature 
above a chosen reference region R.A. 204 58™, Dec. + 39°. The contour brightness temperature unit 
is28°K. The co-ordinates are for epoch 1950.0. The dotted rectangles, A, B, C, D and E, are 
the regions selected for flux comparison measurements (Section 3). 


Tewperature ("%) 


Or gninase 








Fic. 2.—A typical drift curve through the Cygnus X region at Dec. +43° 45’. 


It is assumed that the temperature of this reference region is representative 
of the background temperature for the whole Cygnus X region. Scans at declina- 
tion +49° above the Cygnus X region show that the temperature at a point on 
the galactic plane is not greater than 20°K above the reference point. It is not 
likely therefore that the background flux contributes more than 10 per cent to 
our measurements of flux density. 
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3. Comparisons of 20 cm and 75 cm surveys.—The total integrated flux density 
from a region extending from R.A. 20" 08™ to 20°47™ and Dec. 37° to 45° is 
8600 x 10-7 w.m.~? (c/s)~! at 408 Mc/s. Within the same limits and above 
the same reference region the flux density was measured from Drake and 
Westerhout’s maps and found to be 6200 x 10~*6 and 7400 x 10-*6 w.m.~? (c/s)~? 
respectively. These flux density values are in good agreement thus supporting 
the previous theories that the emission has a thermal origin in optically thin 
HI! regions. 

The spectral characteristics of the individual features were investigated 
both by integrated flux density measurements of five selected areas of the map 
(Fig. 1 and Table II) and by comparison of the brightness temperatures of the 
sources (Table III). 


Taste II 


Comparison at 1390 Mc/s and 408 Mc/s of the flux density from the § selected regions 
shown dotted in Fig. 1 a 


1390 Mc/s 1390 Mc/s 
408 Mc/s Drake Westerhout 
Region R.A. limits Dec. limits x10 wm-? x 10778 wm? x 107° wm7? 
(c/s)- (c/s)~! (c/s)-? 
° , ° , 
+39 30 to +41 I 1g00 1400 1400 
+40 +43 1800 1800 1700 
+43 +45 700 500 600 
+44 +47 600 pate 500 
+42 30 to +46 1000 kes 1000 


TaBLe III 


List of positions and peak brightness temperature of radio sources measured at 408 Mc/s 
together with spectral indices obtained from comparison with Drake (3) and Westerhout’s 
(4) 1390 Mc/s surveys 


Westerhout RA. Dec. 7) peak 17) peak Tos Spectral Remarks 

No. (1950) (1950) 408 Mc/s 1390 Mc/s T4399 index 

h m 8 *) , “K 
63 20 17 30 45 22 104 8 , From Wester- 
hout (4) 

64 20 21 10 37 28 70 ‘ From Wester- 
hout. Poor 
positional 
agreement. 

20 20 00 40 500 ° y Cygni source. 

+048 +2 

20 30 40 39 170 

20 31 30 43 144 

20 33 50 46 47 64 From Wester- 
hout. 

20 33 50 40 178 

20 38 40 41 250 

20 40 20 39 77 

20 42 40 40 70 § 

20 §3 30 43 108 From Wester- 
hout. 


ons] J 
SN Qu Ww 


° 


The probable error of source position (except for No. 66) is 5 min. of arc. 
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From the equal contributions to the flux at the two frequencies from four 
(B, C, D and E) of the five regions it is seen that most of the sources are optically 
thin H 11 regions. This is supported also by the marked similarity of the contours 
at the two frequencies. The ratios of brightness temperatures (Table III) 
show that 7’, varies as f~* for all the sources with the exception of No. 66, which 
lies in Region A (Fig. 1). 

These measurements show that the spectrum of the Cygnus X region is 
thermal throughout, apart from the strongest source No. 66 (the y Cygni source). 

4. The nature of the y-Cygni source.—From the previous section it is seen 
that the ratio of brightness temperatures for the No. 66 source is 18, i.e. T,af-** 
which cannot be accounted for by a thermal origin. Therefore it seems likely 
that at least part of the flux from this source has a non-thermal origin. 

A model is proposed for this source which consists of a non-thermal source, 
flux density «f~! and an optically thin thermal source with brightness temperature 
af-*. Table IV gives the contributions to the brightness temperature at the two 
frequencies from the thermal and non-thermal components in the y Cygni 
source. 


Taste IV 
1390 Mc/s 408 Mc/s 
Thermal Non-thermal Thermal Non-thermal 
22°K 6°K 250°K 250°K 


If there is a non-thermal point source in the position 6f No. 66 which contri- 
butes 250°K to the brightness temperature, then its flux density at 408 Mc/s 


is 230 x 107% w.m.~* (c/s)-*. This represents a practical minimum value for 
the flux density from the non-thermal component since an assumption of a less 
steep spectrum or finite optical depth in the thermal region would increase this 
flux density. If the non-thermal component is appreciably extended, this value 
once again represents a minimum value for its flux density. 

Interferometer measurements at 158 Mc/s on baselines of 220A and 2100A 
show that there is no point source greater than 1 per cent of Cygnus A at the 
position of » Cygni*. Therefore it seems that there is no source less than 20’ 
of arc with flux density greater than 50 x 10-6 w.m.~? (c/s)~! at 408 Mc/s. The 
size of the y Cygni source measured in this present survey is 1°. It appears 
that the non-thermal source has a diameter between }° and 1°. 

From Table II there is an excess flux density at 408 Mc/s of 500 x 10~*6 
w.m.~* (c/s)~! in Region A which includes y Cygni. If this is due to the non- 
thermal source, then to contribute 250°K to the brightness temperature the 
source size would be about 1°. 

5. Conclusion.—The spectral characteristics of the individual sources of the 
Cygnus X region have been studied at two frequencies 1390 Mc/s and 408 Mc/s. 
The results indicate that the spectrum of most sources is that of optically thin 
Hu regions. Thus the emission measures derived by Drake and Westerhout 
from their 1390 Mc/ssurveys are justified. It has also been shown that the spectrum 
of the strongest source in the region is non-thermal and on the basis of the 
relative contributions assigned in the model to thermal and non-thermal com- 
ponents at 1390 Mc/s, it is necessary to reduce the emission measures for this source 


* Data kindly supplied by Dr H. P. Palmer. 
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from 25000 to 20000cm~* pc. Rishbeth (g) tentatively associates an extended 
non-thermal source in the Orion region with an extensive dark nebula. It is 
interesting to note that the y Cygni source also lies in an unusual structure of 
dark nebulosities (3). 
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THE SPECTRUM OF THE GALACTIC RADIO EMISSION 
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Summary 


Absolute measurements of the cosmic background radiation have been 
made at wave-lengths of 1-7 m and 7-9m. _ Scaled aerials were used in con- 
junction with thermally-calibrated receivers in order to make an accurate 
determination of the spectrum of the emission. Observations by Baldwin 
at a wave-length of 3-7 m have been smoothed to provide a comparison at an 
intermediate frequency. 

The spectral index, x, defined by JT, o A*, was found to be 2-37 +0°04. 
There are no significant departures from this value, measured over solid angles 
of approximately 0-2 steradians, except for regions on the galactic plane 
within 60° of the galactic centre. 

The possibility of detecting the extra-galactic component of the integrated 
emission by means of its spectral index is discussed and upper limits are 
derived. 





1. Introduction.—The spectrum of the cosmic radio emission and its spatial 
variation are of utmost importance in investigations of the theories of origin of 
the radiation. Absolute measurements with large and complex aerial systems 
have proved very difficult. A comparison of such measurements made at 
different frequencies by various workers yields a wide range of values for the 
spectral index. 

A more promising approach is to use identical aerial systems at various 
frequencies (1). Direct comparison of aerial temperatures may then be made 
and the problem is reduced to a determination of receiver sensitivity. The 
use of scaled aerials necessarily limits the resolving power that can be easily 
obtained. It does, however, permit comparison of the more extended com- 
ponents of the galactic emission, and the mean temperatures derived may be used 
to calibrate larger instruments. 

This paper contains the results of measurements of the brightness of the 
northern sky at wave-lengths of 1-7 and 7-9m. Corner reflector aerials with 
reception patterns of 15° x 40° were used. 

The magnitude of an isotropic extra-galactic component of the background 
has been discussed in relation to a number of cosmological models (2, 3). 
Observationally, such a component is indistinguishable from the outer shell of 
the galactic halo if it has the same spectrum; an upper limit may be placed on 
its magnitude from considerations of symmetry (4). If, however, its spectrum 
differs from that of the galactic emission, the present observations allow somewhat 
finer limits to be set than those determined above. 





No. 3, 1960 The spectrum of the galactic radio emission 249 


2. Apparatus.—The aerials used in these observations were of the go° corner 
reflector type with sides 1-1 by 4A. Each contained 4 full-wave dipoles 0-45A 
from the apex. ‘They were mounted on an East-West line, 0-2A above the ground 
and could be tipped in declination. The reception pattern of each aerial covered 
15° in right ascension and 40° in declination. They were matched prior to the 
observations to give a S.W.R. < 1-05. 

Phase-switching receivers (§) were used at both frequencies. A ‘ hybrid-ring’ 
network provided the two correlated inputs required for this type of receiver. 
It is shown in the next section that this is effectively a Dicke (6) system but has 
important advantages when comparing noise sources of differing impedances. 

Two sources of noise power were used to measure the sensitivity of the 
equipment: a temperature-limited diode source for daily measurements, and a 
thermal noise source for the absolute calibrations. The thermal sources used 
have been described elsewhere (1). 

3. Observational method.—The arrangement of the apparatus is shown in 
Fig. 1. 












































Fic. 1.—A schematic diagram of the apparatus. 


The aerial or noise source is connected at A through the matcher to B. The 
comparison source at D is a thermal load immersed in water in a thermos flask 
and is maintained at room temperature. ‘The matcher is adjusted so that the 
impedances presented to B and D on the network are equal. 

Let the noise powers at B and D be represented by their effective temperatures 
T,, and T,,, and that of the receiver by T,. The temperature at the input of the 
receiver in the ‘in-phase’ position of the switch is given by 


i 
T+ (=+ 2) +(- =) T,+T ° (1) 
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The corresponding expression for the anti-phase position of the switch is 


T,+ ( + =") < (3 - =*) =Ty+Tp. (2) 

The output of the phase-switching receiver is proportional to the difference 
between (1) and (2), i.e. is proportional to 

T, ne Ty. 

The system is therefore equivalent to the more straightforward Dicke system 
where the receiver is switched to the aerial and the comparison source in turn. 
Its advantage lies in the fact that a C.W. signal (at the centre frequency of the 
system) introduced at E, appears at C attenuated by an amount proportional to 
Z,—Zyp, where Z, and Z, are the impedances presented to the network by the 
matcher and thermal load respectively. With a carefully constructed network, 
the matcher may be adjusted to give an insertion loss between C and E of the order 
of 80db. The matched condition, when the various sources are interchanged, 
may therefore be reproduced with great precision. 

The arguments above refer to a single discrete frequency, whereas the obser- 
vations are made of noise signals with systems of finite bandwidth. _ It is essential 
that the matching be done at the effective centre frequency of the noise band, 
otherwise the receiver output may be dependent on the impedance of the noise 
source. This condition may be checked by transforming the source impedance 
through an arbitrary length of cable. After re-matching, the receiver output 
should be the same as before. 

The procedure adopted may be summarized as follows : 

(i) A large modulated C.W. signal is introduced at E. 

(ii) ‘The matcher is adjusted with the aid of the audio-detector to reject the 
modulated signal. 

(iii) The receiver is reconnected to E and the recorder output noted. 

(iv) Steps (i) to (iii) are repeated with an additional length of cable (usually 
A/4) between the noise source and the matcher. 

(v) If necessary, the frequency of the C.W. signal generator is adjusted and 
the procedure repeated until introducing an arbitrary length of cable and rematch- 
ing produces no change in the receiver output. 

The aerials were moved in 20° steps between declinations +20° and +80’. 
At each setting a continuous record was taken over a period of 24 hours. Each 
run was repeated until at least two records, consistent throughout to better than 
5 per cent, were obtained. Daily measurements of sensitivity and zero were 
made, some through the aerial cable in order to determine the cable attenuation. 
A comparison of the diode noise source with a thermal load had been made by 
F. G. Smith prior to this set of observations. This comparison was repeated by 
the author. 

4. Method of analysis.—The records selected as described above were sampled 
at hourly or half-hourly intervals. The measurements of cable attenuation, 
zero-level, and sensitivity were used to convert the record deflections to equivalent 
aerial temperatures. Where a number of values were available, the mean was 
computed. The results are shown in Figs. 2, 3, 4 and 5. 

The results of a survey by Baldwin (7) at 3-7m with a system of higher 
resolving power were smoothed over 15° x 40° in order to permit a comparison 
with the present observations. ‘The smoothed scans are also shown in Figs. 2-5. 
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Fic. 2.—The variation of aerial temperature on Dec. +20° at wave-lengths of 1°7 m (continuous 
line), 3°7 m (crosses) and 7°9 m (filled circles). 














Fic. 3.—The variation of aerial temperature on Dec. +- 40° at wave-lengths of 17 m (continuous 
line), 3°7 m (crosses) and 7-9 m (filled circles). 


The spectral index was then determined by finding the scaling factor which 
gave the best agreement between the temperature curves obtained at the three 
frequencies. No further adjustments of temperature scales have been made. 
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The effects of aerial smoothing and reception through subsidiary maxima may 
be neglected when comparing the 1-7 m and 7-9 m observations. The smoothing 
function used for the 3-7 m observations contained no side lobes. The highest 
temperatures are therefore too high and the lowest very slightly too low in relation 
totheothercurves. This effect may be as much as 5-10 per cent near R.A. 19" 00™ 
on the curve for Dec. + 20°, but it will be much smaller for the other parts of the 
sky. 














Fic. 4.—The variation of aerial temperature on Dec. +-60° at wave-lengths of 1-7 m (continuous 
line), 3:7 m (crosses), and 7-9 m( filled circles). 














Fic. 5.—The variation of aerial temperature on Dec. +80° at wave-lengths of 1-7 m (continuous 
line), and 3-7 m (crosses), and 7-9 m( filled circles). 


The other sources of error that must be considered are listed below: 

(i) The loss in the cable between the aerial and the receiver. ‘This quantity 
was measured a number of times and a mean value obtained with an accuracy of 
I per cent. 

(ii) Variation in the temperature of the cable with respect to the reference 
thermal load. With the cable attenuations concerned this effect would be of the 
order of 1 per cent and has been neglected. 

(iii) The reading accuracy of the records. This is about 1 per cent for one 
individual record and the error of the mean is therefore insignificant. 
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(iv) Fluctuations in receiver sensitivity and zero level. The receiver sensitivity 
was maintained at a nearly constant level by means of an AGC system. The 
7°9m observations were often seriously affected by daytime interference from 
R.A. 17"00™. As mentioned above, records showing departures from the 
mean curve of more than 5 per cent were discarded so the relative accuracy of 
most of the scans should be about 2-3 per cent. Where bad observing conditions 
persisted (00" to 06", Dec.+ 40° at 1-°7m; 17" 00™ to 24" 00”, all declinations at 
7-9 m) wider tolerances had to be accepted and the accuracy is of the order of 5 per 
cent. ‘This is also true for Dec.+80° where only 2 scans were made at each 
wave-length. 

(v) The accuracy of the noise power standard. ‘The calibrations of the diode 
noise source by the author agreed with those of F. G. Smith, to within 2 per cent. 
An absolute accuracy of 5 per cent is assumed. 

It is estimated that the temperature ratios obtained from most of the individual 
points are accurate to about 10 per cent. Where a large number of points are 
used, the accuracy improves to about 7 per cent, which is the present estimate 
of the absolute accuracy of the whole system. 


5. Discussion. . 


(i) The spectrum of the background radiation.—The most obvious feature of 
the curves shown in Figs. 2-5 is the excellent agreement at all three frequencies 
over most of the sky. It is clear that the components responsible for a very high 
proportion of the integrated emission all have the same spectral variation. Only 
on the galactic plane at Dec. + 20°, where the effects of H 11 clouds may be expected 
to become important, is there any clear evidence for variation in the spectral 


index. ‘The observations therefore confirm the conclusions of Adgie and Smith 
(x) based on observations over the same frequency range with aerials of lower 
resolving power. 

Another feature of probable significance is the lower intensity at 1°7m from 
R.A. 12"co™-16"30™ on Dec. +20°. In this region, a substantial portion of 
the radiation comes from a narrow belt of emission extending from the galactic 
plane near /=o, up towards the galactic pole. A similar reduction of intensity 
at 175 Mc/s is also shown in the observations of Adgie and Smith (1). Additional 
support for a steeper spectrum in this region is provided by the survey of Kraus 
and Ko (8), where the feature is much less prominent than at lower frequencies 
(7, 9, 10). 

The absence of a significant variation in the spectral index over most of the 
sky allows the observational errors to be averaged out. The accuracy of the 
spectral index is then determined primarily by the absolute intensity scales. All 
points oo" oo™ to 17° 00™ in Figs. 2—5 have been replotted in Fig. 6. The value 
of the spectral index determined from these points is 

X¥ = 2°37 + 0°04. 

This value is also consistent with the 3-7 m observations. 

(ii) The extra-galactic component.—The determination of the integrated 
emission from extra-galactic sources and hence the average emission/unit volume 
of extra-galactic space is of considerable interest, both as a possible means of dis- 
tinguishing between different cosmological models (2, 3) and as evidence on the 
nature of radio stars (11). The present observations now make it possible to 
obtain additional information. 
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Part of the extra-galactic emission is due to normal galaxies, whose spectra 
may be expected to be similar to our own; observations of the Andromeda 
nebula (12, 13) have given figures close to that found in the present observations. 
Part will however be due to the more powerful extra-galactic sources; all of those 
whose spectra have been measured have been found to have an emission which 
varies as A**—A!?, The integrated emission will therefore have a component 
whose brightness temperature will lie in the range A**—A**. An estimate of the 
upper limit of this latter component is of great interest (11). 

Whilst the incidence of radiation from normal galaxies will appear as an equal 
displacement at all wave-lengths of the zero level of the galactic emission (the 
dotted line in Figs. 2—5 indicates the upper limit of the total extra-galactic emission 
derived by Baldwin (4)), extra-galactic emission having a steeper spectrum will 
give rise to unequal displacements at different wave-lengths. 


178 Mc/s 





L 


8 
Ty (°K 210°) 





Fic. 6.—The aerial temperature at 1-7 m plotted against the aerial temperature at 7-9 m. The 
solid dots are from Dec. +20°, R.A. 13%00-16400. 


This effect is shown more clearly in Fig. 6. If an extra-galactic component 
with a steeper spectrum were present, the line through the points would no longer 
pass through the origin but would intercept the 7-9 m axis. In fact, a line fitted 
to the present observations by the method of ‘least squares’ does intercept the 
7°9m axis at 7.,=400°K. In view of the possible errors in the absolute intensity 
scales, this value cannot be regarded as significant. A maximum value of goo °K 
at 7°9m has been determined for the excess emission. Table I shows the maxi- 
mum possible magnitude of any component of the integrated extra-galactic 
radiation due to sources with a spectral index in the range S oc A°* to S oc Al®, 
corresponding to brightness temperatures proportional to A** to A*?. 


TABLE I 


Extra-galactic contribution Spectral index 
(A= 3-7 m) 
<240 “K 28 
<155 K 30 
<1ie °K 3°2 
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The contribution from normal galaxies has not been included. The last 
value corresponds to the mean spectral index determined by Whitfield (12) for 
63 radio sources at galactic latitudes greater than 10°. 

I would like to express my thanks to Dr F. G. Smith, who was largely respon- 
sible for the development of the experimental technique and who also assisted 
with the 7-9m observations, and to Dr J. Baldwin who assisted with the 1-9m 
observations. My thanks are due also to Professor M. Ryle for helpful discussion. I 
am indebted to the National Research Council of Canada for a Special Scholarship. 

Mullard Radio Astronomy Observatory, 

Cavendish Laboratory, 
Cambridge : 
1959 June 19 
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A COVARIANT FORMULATION OF THE LAW OF 
CREATION OF MATTER 


F. Hoyle 
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Summary 


A covariant law is given for the creation field of matter. When the 
electromagnetic and nuclear fields, and the stress within matter, are omitted 
from the energy-momentum tensor, the Einstein equations possess an infinity 
of isotropic, homogeneous, stable, steady-state solutions. The possibility 
exists that, if the neglected terms were reinstated, the solutions would possess 
a slow secular drift. 


1. Introduction.—The object of the present paper is to formulate the steady- 
state model of the universe in a satisfactory covariant way. A brief review will 
first be given of the existing state of affairs, so far as theories based on the Einstein 
field equations are concerned. For considerations lying outside this framework 
the reader should consult Bondi and Gold (1). 

Isotropic, homogeneous models of the universe possess a line-element of the 
form 


r 


ds? =c2 dt? RX) ar" 4. 19(d0* + sin® Odd? )| (1) 
1 —kr? 
where k can be zero or +1. The well-known first-order approximation to the 
red-shift magnitude relation follows already from (1), and therefore does not 
depend at all on the system of dynamics used to determine the precise variation of 
R with time ¢. 
The dynamics of cosmological theory are determined by Einstein’s equations, 
o, . hg. G + Ag, =-—K ae (2) 
where G,,, is the Ricci tensor, 7',, is the energy-momentum tensor, x = 87ry ce! with 
y the Newtonian constant, and A is the cosmological constant. For a single-stream 
fluid with isotropic pressure, 7" contains the term 


(p+p)o*v" — pyr’, (3) 
where vo” is the velocity dx*/ds, p is the material pressure, and p is the density 
measured by an observer moving with the fluid. (Even for this observer, the 
particles of the fluid are in random motion, so that the apparent masses of the 
particles differ from the rest masses. ‘This ditference must be included in p.) 

Expression (3) does not comprise the whole energy-momentum tensor. 
Extra terms containing the electromagnetic field and the nuclear field must also 
be included. But these are usually taken to be negligible, and (1), (2), together 
with (3) for 7’, define the well-known forms of relativistic cosmology described 
very clearly and in detail by Robertson (2). 


‘ 
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A theory containing the creation of matter can be obtained (3) either by 
introducing a suitable term into (2) or by the entirely equivalent procedure of 
adding an extra term to 7*’. As was pointed out by McCrea (4) the second way 
of writing the theory is preferable because it emphasises the crucial point that the 
logical structure of the Einstein equations is not thereby changed. 

Equation (2) connects a geometrical tensor of defined form on the left with an 
undefined physical tensor on the right—the structure of the relativity theory does 
not itself decide what form the physical tensor shall take. (In fact the formulation 
of T+’ has changed since Einstein's paper of 1915, in the sense that the nuclear 
field must now be included in 7“’.) In particular, conservation equations, 
T*’, v=o, hold whatever the formulation of T+’, . 

For a universe of infinite volume (as in the steady-state theory) energy con- 
servation for the whole universe is a meaningless notion. Conservatior of energy 
must be considered in relation to a box of finite volume. ‘Iwo important cases 
evidently arise: the box can have a fixed proper volume, or its walls can expand as 
the universe expands. When there is creation of matter, the conservation equa- 
tions require energy to be constant in a box of fixed proper volume. When there 
is no creation of matter, the conservation equations require energy to be constant 
in a box that expands with the universe. 

To end these introductory remarks, it may be noted that if the time-axis is 
open to — ©, as it is in the steady-state theory, the cases k=o, +1 are indistin- 
guishable over any finite portion of the universe. Hence for any comparison with 
observation it is sufficient to consider the case where & in (1) is taken as zero. The 
theory (3) was further simplified by putting A=o. ‘These simplifications will be 
retained in the following work. 

2. The creation term.—'Yhe success of the isotropic, homogeneous theory in 
predicting the linear red-shift magnitude relation is a strong indication that the 
form (1) of the line element is correct (apart of course from detailed fluctuations 
due to local irregularities in the distribution of matter). Moreover, observation 
provides direct evidence that the large scale distribution of galaxies is in fact 
isotropic and homogeneous. 

This raises the remarkable point that apparently a preferred direction is defined 
at every point of space-time, namely the direction of the universal time ¢, viz. 
(1, 0, 0, ©) in the system of coordinates used in (1). By differentiating this 
direction covariantly, a tensor C,,, is obtained, of which 

Cio = Ca, = Cag = — RR (4) 
are the only non-vanishing components. In essence, the procedure adopted in 
(3) was to multiply C,,, by a constant « and to add aC’ to (3) in order to form 7%’. 
The resulting dynamical equations have a steady-state solution 

R(t) = exp (ct/a), 
p= 307/8:ya’, (5) 
a= 30 /8ray. 
It was also shown that any fluctuation of R(t) from this solution is exponentially 
damped in a time interval of order a/c. 

In this way a unit of distance a, a mass density p, and a time unit a/c, were 
obtained. ‘The emergence of universal units of length, mass, and time arises from 
the introduction of the creation constant «, whici, .ogether with c, y determine a 
andp. The time unit a/c is the reciprocal of the Hubble constant H. 

18 
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The objection to this theory—that the definition of C,, makes appeal to a 
particular coordinate system—was answered (5) with the argument that no theory 
working entirely in terms of covariant laws (i.e. without reference to a special 
coordinate system) could explain the development of preferred directions in space- 
time. Underlying this argument was the concept that the ultimate aim of cosmo- 
logical theory should be to establish the following result: given any arbitrary 
distribution of matter and motions, and arbitrary values of the metric tensor and 
its derivatives, on a space-like surface, prove that the universe must ultimately 
evolve to its present isotropic, homogeneous state. It is difficult to see how this 
result could ever follow from covariant laws, since neither the laws nor the initial 
situation single out any particular coordinate system for especial preference. 

The point is so important as to be worth restating in a somewhat different way. 
One learns from the relativity theory that space and time are no more than a 
fourfold set of parameters by which physical events are ordered and that the 
parameters can be chosen in an infinity of equivalent ways. Experience in micro- 
scopic physics provides overwhelming evidence in favour of this point of view. 
Yet special coordinate systems arise in the macroscopic world, not only as a matter 
of convenience, but in a quite fundamental way. Newtonian theory is dependent 
on the concept of absolute space. Ata more sophisticated level, Mach avoided the 
necessity for absolute space by choosing a coordinate system that was non-rotating 
and unaccelerated with respect to the distant bodies of the universe. At a still 
more sophisticated level, in the theory of relativity the behaviour of an isolated 
singularity in accordance with the Schwarzschild solution is dependent on the 
existence of a suitable isotropic boundary condition at large distance from the 
singularity. It is of course the isotropic character of the universe that supplies 
this boundary condition. Indeed it is in just this way that the special coordinate 
system determined by the large scale structure of the universe impresses itself on 
everyday experience—for example, in controlling the development of cyclones 
in the atmosphere of the Earth. 

The dilemma that exists between macroscopic experience on the one hand and 
the covariant outlook of microscopic physics on the other hand is therefore very 
real and serious. ‘The dilemma is side-stepped by an unsatisfactory device in 
systems of cosmology where the negative time-axis is closed. In such cosmologies 
the universe is supposed to have been created at a finite value of t. In the mathe- 
matical sense, ‘ creation ’ means that certain boundary conditions are imposed 
on the differential equations of physics. Isotropy and homogeneity simply make 
their appearance at the outset among these boundary conditions. 

When the time-axis is open to — 00 the question cannot be dismissed so easily — 
or so artificially. ‘Two points of view can be put forward. One of these has 
already been mentioned—of starting with an entirely arbitrary distribution of 
matter and of proving the development of isotropy and homogeneity. 
This point of view apparently demands the adoption of some non-covariant 
law. 

Alternatively, one can argue that such a demonstration would be quite un- 
necessarily strong, since all that really needs proving is the stability of a large scale 
isotropic homogeneous distribution. According to this second point of view, 
the notion of starting with an arbitrary distribution of matter could well be an 
invalid concept, because the universe need never have been in any state other than 
one of homogeneity and isotropy however far back in time we go. 
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Although the choice of procedure is to some extent a question of philosophy, 
the possibility that the law of creation can be made entirely covariant in the second 
case is evidently of physical significance. ‘The problem therefore arises in the 
following form: how shall a covariant law of creation be formulated that maintains 
permanent large-scale homogeneity and isotropy? In the following section a law 
will be obtained that gives timewise stability. Spacewise stability is not strictly 
proven, but the very form of the law makes spatial stability extremely probable. 

3. A covariant creation law.—It is proposed that matter be thought of as 
radiating a creation field. For simplicity, the field will be takem scalar, thus 

C|¢ as Kp, (6) 
where K is a constant, p is the mass density, and [] implies that the field @ is 
differentiated twice covariantly and that the two suffices are then contracted. This 
tield presumably arises, if it exists at all, from the microscopic pfocesses of funda- 
mental physics. Although at present nothing is known of what these processes 
might be, it is reasonable to suppose that the source of the field would be propor- 
tional to the mass density p. Uncertainty concerning the microscopic origin of the 
field is therefore limited to uncertainty in the value of K. (The present formalism 
includes the case where there is no field ¢, this being given by K =o.) 

Differentiation of g gives a preferred direction at every point where the 
gradient is non-vanishing. A preferred direction could also have been obtained 
from a vector field, but there would be a far greater uncertainty in the nature of 
the source function in the case of a vector field. ‘Thus for a field gp, the equation 
corresponding to (6) would be 

9, =s K, p, 
and a considerable knowledge of the microscopic processes would be needed to 
determine the direction of K,, at each source particle. 

To obtain a system of cosmological dynamics, C,, is first obtained by differenti- 
ating @ twice, viz. 

" o* 0 
Cu se Tb @) 
where I’,) is the usual Christoffel symbol. ‘The tensor C,, is then added to (3) in 
order to form 7,,,. (On this occasion nothing is gained by multiplying by a 
constant before C,,, is added to 7, ,, since any such constant can be absorbed into ¢. ) 

It is seen from (7) that C,,, is of necessity symmetric in its two suffices, so 
that the symmetry of 7”, is not impaired by this procedure. 

Consider now the dynamical equations for the smoothed case in which @ is a 
function of t only. With k=o, the line element (1) can be written in terms of 
coordinates f, x,, X», X3, very simply as 


ds? = c2 dt® — R®(t) (dx,2 + dx,? + dx;?). (8) 


The only non-vanishing Christoffel symbols are 
n° = RR5,,/c?, 15; == R8,,/R; i,j=rt, 2, 3, 
while the only non-vanishing components of G,,, are 
Giy= —(RR+2R)8j/2;  — i,j=1, 2, 3, 
Goo = 3R/R. 
From (10) it follows easily that 


G = 6(RR + R*)/ R®e?. 
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The non-vanishing components of C,, are 
Cy=—RRpbyic; i, j=1, 2, 3, (12) 
Coo = y- 
Three independent dynamical equations follow from (2), (7), and from the 
construction of 7,,,. With p =o in (3) and taking v = c(1, 0, 0, 0), these equations 
are 


2RR+ R*—KRR¢=0, (13) 
(14) 


39.5 = Kep—$. (15) 

Consider now the condition for a steady-state solution, R= exp (Ht) with Ha 
constant. From (13) we require 

G=3H|k, ie. p=3Ht/K, «K=B8ry/c. (16) 
Thus ¢=0 and (14) gives 
p= 3H?/8my, (17) 
as in the previous theory (3). Substituting for ¢, p in (15) now gives K = 3c*. 
‘Thus unless (6) takes the form 
C]p = 3e%p, (18) 
a steady-state solution cannot be found at all. 

This result is most satisfactory since it removes the former arbitrariness of a 
multiplying constant (the constant « in section 2). ‘Thus the steady-state theory 
determines the precise law of the creation field, viz (18), thereby making a definite 
statement as to the form of the law that microscopic physics should reveal. 

The disappearance of the multiplying constant « means that universal units of 
length, mass, and time are no longer defined, however. The difference shows 
itself in the fact that, provided @ satisfies (18), a steady-state solution of (13), (14), 
and (15) can be found for any value of the Hubble constant H. A discussion of 
whether this is an advantage or not will be deferred until after the timewise stability 
of the equation has been discussed. 

It may be observed that the form of (18) requires g to be determined always 
by the distant matter of the universe, so that ¢ is not appreciably affected by local 
fluctuations in the distribution of matter. For this reason it is to be expected that 
the universe is spacewise stable. 

4. Stability with respect to time.—-Suppose (18) is satisfied. Consider a per- 
turbation from the steady-state solution R=exp(H,t), where H, is a constant. 
‘The perturbation can be expressed by 

: Hy, (1 +€), €€1. (19) 
Differentiating (19) gives 
2 
: 7 = ? Hyé ° (20) 
By substitution for R, R in (13) we then have 
xb = 3H, (1 + €) + 26, (21) 
asecond order term in 3€ being neglected. Putting K = 3c* in (15), and eliminating 
p between (14) and (15), gives 


: + 4x. (22) 
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By differentiating (21) to give 
Kp = 3H, € + 2é. 
and by substitution for R/R, @, g in (22), we finally obtain 
42+ 9H f=0, 
with a general solution of the form 
£= A exp (—gHot/4) + B. (24) 

Hence we see that a small perturbation of the universe introduces a damped 
term, together with a constant shift of the Hubble constant. ‘The possibility of 
such a shift is obviously to be expected, since a steady-state solution can be found 
at any value of H. ‘The interest therefore lies in the stability of the time dependent 
term. ‘The universe evidently has no exponentially unstable solution. Nor is 
there any solution in small oscillations. ° 

5. Discussion.—The point of view of the present theory is that the energy 
momentum tensor 7\,, contains four parts, viz. 

Tee Tet Tet T+ T,”", (25) 
and that (25), together with the Einstein relation (1), determines the dynamical 
equations of cosmology. ‘The quantity 7',"", given by (3) for a single stream fluid 
with isotropic pressure, generates the ordinary gravitational field. ‘The electro- 
magnetic energy-momentum tensor 7',“" is given by 

To" = 1/40 (— Fe Fy + ge F’ F,,), (26) 
where F" is the electromagnetic field tensor. ‘The energy momentum tensor of 
the creation field is given by 

lp = 3¢%p, 

T=", (27) 
the covariant form of 4” being written out explicitly in (7). No canonical form 
for the nuclear energy-momentum tensor 7',”” is at present known. 

In this way the creation field plays an entirely similar role to the other fields. 
From a cosmological point of view the important terms in (25) are 7,*”, 7,"". 
Ironically, it is 7", 7|,“" that are studied so extensively in microscopic physics. 

It is possible to take two views of the present value of the expansion rate R/R. 
Without perturbations, and accepting (27), the universe is steady, in which case 
there need never have been any other value of the expansion rate. It is 
then meaningless to attempt to find a reason for the observed value, 
R R= 2:5 x 107*sec™!. 

On the other hand it is possible that perturbations, satisfying (24) with Bo, 
arise from 7',“", 7.,“", or simply from local spatial fluctuations. In such a case the 
universe would possess a slow secular drift in its expansion rate H. Consequently 
the mean density of matter, given by (17), would also possess a secular drift. 

Whether or not such a drift would actually occur is excessively difficult 
to determine. ‘The main contribution to 7,“” comes from starlight. ‘The main 
contribution to 7',“” presumably comes from nuclear binding and from the neutrino 
field (these have been lumped together into one term). It should be a good 
approximation to regard these tensors as spatially smooth, in which case it is to be 
expected that the constant K in (6) can be chosen so as to permit a stable, steady 
solution R= exp (Ht) of the dynamical equations, provided the tensors T,*", T,,”" 
are steady. ‘This is the doubtful point, however. ‘The formation of galaxies, and 
the incidence of nuclear activity inside stars, must be affected by the value of H. 
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Hence a secular drift of H must change T,"’, 7,,“’, which in turn produces feed- 
back into the dynamical equations themselves. If the feed-back were positive, 
drift would continue. 

In so far as any criterion can be given for judging a system of cosmological 
dynamics—other than the overriding observational criterion—it might be this: 
that the ‘ best ’ system of cosmology is the one that admits of the greatest variety of 
physical possibilities. If this criterion is admitted, a universe with slow drift is 
to be preferred to an entirely steady-state universe at fixed H. In the latter case 
only one unique set of properties can be worked through, namely those that 
correspond to the one unique value of H. In the case of secular drift, on the other 
hand, an infinity of different sets of properties can be examined; for if the drift 
is slow enough the full physical implication of each value of H can be established in 
turn. (Since different H values imply different values of the mean density p, any 
marked shift of H certainly implies a non-trivial change in the properties of the 
universe.) Defining H~! as the length of a ‘ generation ’, with H always taking 
its instantaneous value, it is possible if the drift were sufficiently slow to proceed 
backwards along the time-axis through an infinity of generations and yet for the 
universe never exactly to have repeated itself in any two of these generations. 


St John’s College, 
Cambridge : 
1959 August. 
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Summary 

It is shown that the Hund gravitational equations can be applied to a 

spherically expanding universe. The experimentally determined constants 
are used to measure the age of the universe. 





1. Introduction.—A variety of cosmologies have been invented to describe 
the nature of the universe. Any cosmology must, in order to be considered as 
representing the universe, be consistent with two macroscopic properties of the 
physical universe. First, the mean density of the universe is constant over space 
within measurable distances and, secondly, the velocity of recession of nebulae 
(the red shift) is proportional to the first approximation to the distance from the 
observer. The object of this paper is to provide a semi- Newtonian cosmology 
for the universe, and from this to determine the age. 

Hund (1) has suggested that gravitational fields may be expressed in terms 
of a polar vector f and an axial vector g which act upon any gravitational mass m 
with the force 

m{f + (v/c) x g} (1.1) 
v being the velocity of the mass. The fields f and g satisfy the following differ- 
ential equations. 


V-g=0 (1.2) 

(1.3) 

V - f= —4ryp t (f* +8*)/(2c?) (1.4) 
(1.5) 


p is the mass density. 

The fields f and g are supposed to be produced by gravitating masses and 
currents, analogous to the electromagnetic field. All forms of energy and energy 
flux, including the fields due to f and g themselves, are considered as sources of 
gravitation and so those terms in the field equations which correspond to the 
source terms in Maxwell’s equations must include in addition terms involving 
f and g. 

The rate of creation of mass per unit volume is derived from these equations 
to be 

ep 
at 
and this is the rate at which gravitation does work on the system divided by c?. 


+9. (pv)= =pf.v (1.6) 
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These equations may be used to provide a smoothed out representation of the 
universe. If it be assumed that everything is radially symmetric g vanishes and 
f and v are in the radial direction. 


Under these conditions equations (1.2)-(1.5) reduce to 


_ 


Cc 


£2 f= —4mp+ frllae) (1.7) 


or 
Equation (1.6) becomes 
3 


5 = (r*pv)= “hep. (1.9) 


~ 
Co 


+ 

ct 
Furthermore, if it be assumed that the motion of the universe is generated solely 
by its own gravitational interactions, f will be the acceleration of the matter under 
consideration and will be related to v by 


(1.10) 


d : — ; 
where ai denotes differentiation following the motion. 


Equations (1.7), (1.8) and (1.10) are three differential equations relating 
the three quantities f, v, p. 

Having solved these equations, the quantity (fvp)/c?, the rate of creation of 
mass per unit volume may be derived. It follows that in this theory there is 
bound to be creation as (fvp)/c? cannot vanish except for the trivial zero solutions 
of equations (1.7) (1.8) and (1.10). This, as mentioned previously, is equivalent 
to the rate at which gravitation does work. The actual solution of equations 
(1.7) (1.8) and (1.10) is clearly an almost impossible task. It is however possible 
to find approximate solutions which hold within certain ranges. 

2. Application to a class of model universes.—If p is uniform, Hund (1) showed 
that the solution of (1.7) which is finite at r=o0, is 
\ 27ypr) 


3 =} 


facts 
lr 
Now if (V 27yp r)/c is small 
f= — 4rypr. (2.2) 
This expression will be used to determine an expression for the time variation of 
density. This will be legitimate provided (4/27yp r)/c is small. 
Using equation (1.8) 

1 dp _ 
: % 
where 7+ is the reciprocal of Hubble’s constant (2). It will be seen that wv is 
proportional to r._ Furthermore because there is recession v is positive and 


dp is negative 
at gi . 


rir at present (2.3) 


Substituting equations (2.2) and (2.3) in equation (1.10) 


(4) ae a4. } I dp 
p\at}j | * padtf\ zp dtl’ 
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This simplifies to 


This may be rewritten as 


dp\*\ _ -2'3 
+) } = 8:ryp*, 


2 
(F) = 24myp* + Kp*® 


If py is the present density, it follows that 


I =) 8 fp ( 3 ) (2) 
— — = —TYPy <— Ss a — J — 
(5 dt 3 ” (Po Sarypyt* Po f 


and the nature of the solution depends on the sign of 


This integrates to 


3/(82rypot?) —1 =A say. 


Writing (p/p))'* =x, equation (2.7) becomes 


(: ai) = 87yPq {x3 + Ax), 
x dt 3 
The behaviour of x is qualitatively similar to that of p. 
(a) If A is zero the solution of equation (2.4) is 
I Poly” 
PO” bry(t+h) (t+t)? (2-9) 
if t=0 be the present epoch. 

This implies that the density was originally infinite at a time ty = (67yp,)~!? 
back and will decrease steadily becoming zero at infinite time. f, is thus the age 
of this model universe. 

(b) If A is positive (=a say) equation (2.8) can be rewritten as 


i. | a en 
3 F | 3070" 


‘a+x—yva +i—v/a 
SANE) ~ log (18) 
Vatrx+rvya Yari+vya 





Vat Vari 
ax a : 





Here again x (and hence p) decreases steadily with time. 
At time t, past, where 


8 2 a+l I a+i-a 
(=*") gw MAS? xalog (+221 —¥*), (2.12) 
3 a 2a” Va+i-va 


x (and p) were infinite. If a tends to infinity, ¢, tends to r. ft, is thus the age of 
this model universe. At infinite time in the future x (and p) will be zero. 7 is 
thus an upper bound to the age of the universe. 
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(c) If Avis negative (A b say). 


Siryp,\'? [ dy : 
of bas {= ee I>xy>h ‘ a. 
( 3 ) J) yy -5) (y>x2b) (2.13) 


The situation here is slightly different. 


x (and p) decrease steadily with time 
until x=bwhent=?'. 


After that time x (and p) will increase again until x becomes 
infinite. ‘This takes a further time t’ +t, where ¢, is the time past at which x 
was originally infinite and is thus the age of this model universe. 


eS) sd , ( dy 
3 Sa (yy —5)! 


. 


hyi2 
a cos pra} ( he (2.14) 


The time f' is given by 
[ dy 
Jo y*(v—6)'? 


hyi2 
(1 — 753008 1 pie 
, iad 


+ 
b) 


12 
Hence $. + bd : 2.16 
ence (fy +0) m2] (2.16) 


8 Ty 12 
re ty. (2.15) 


is thus the lifetime of the universe from its inception at infinite density, through 
its expansion to a density p,b*, to its cessation at infinite density again. 

It may be remarked that this analysis holds only provided that (V 2zypr/c) 
is small, and as p becomes infinite, the corresponding value of r becomes zero. 
Thus the solution is not strictly true at the limits —t, and t’, but will give an 
idea of the manner in which these limits vary with the present density of the 
universe. 

3. Application of observational values.—The presently accepted value of 
Hubble’s constant is given by (2). 


7-l=2:5x107sec™!, += 13 x 10° years (3.1) 
and it is estimated that the present mean density of the universe is between 
10°" and 1o-*'gm/cc (3). (3.2) 
Now if the present mean density is 1:08 x 10~** gm/cc 
3/(8mypor*) = 1 (3-3) 


and A vanishes, and so in order to discuss the effect of mean density, the following 
values of p, will be taken, which give a fair spread. 


Taare | 
Po (a@m/ec) 3/(Srryper*) A 
1°08 x 10~° 10 ° 
1o8 X10 ' ‘a 
r'o8 x 107° “T 
1'0o8 X 10 10-4 


== o’9 
— 0°9999 
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Case (a): py= 1°08 x 10° gm/cc 


Srryp,\! 2 V 10 I anes 
petal i an. 28 «ne _— l r7(\ o-+ 9 . 
( 3 ) 0 Ps 25 og(v 1 3) (3-4) 


which may be rewritten as 


V10 

t,/t = hie log (V 10+ 3) (3.5) 
=e 
whence 
ty = 0°'90T (3.6) 
= 11°7 x 10° years. (3.7) 
Case (b): py = 1:08 x 10>** gm/ce 

ty = 1/(67rypq)'* (3.8) 
ir (3-9) 
8-7 x 10” years. (3.10) 


Case (c): pp = 1°08 x 10°* gm/ec 
Po £ 


Smryp,\"2 Saryp,\"2 
(“™20) te ( TY f *) de. at 
3 a tite Ok 


where 


(men (1-0) + nein 4). (3-12) 


3 h p3/2- 


Now it will be seen that as we are only interested in orders of magnitude we may 
assume (1—6)!? small. (If p, be one degree of magnitude smaller, the formulae 
(2.13-2.16) do not apply, whereas if p, be greater the approximation is a fortiori 
better. ) 

Taking then terms to the first order only in (1 —6)'?, (3.12) becomes 


- 12 
(=) 1’ =2(1—b)' 
i 3 
12 
. ( pe te 
mz) 


t'/7 = 3/(47ypoT*) 


ty 7 3 is 6’ 
tT 2\8rypyr, r 


For the value of p,) under consideration 3/(87ypy)r*) = 0-1 and so 


and 


and from (3.11) 


t'=0°20T 
= 2°6 x 10° years 
ty = 0°30T 


= 3°9 x 10° years. 
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Case (d): py = 1-08 x 10-” gm/cc. 
In this case 
3/(82rypor*) = 10-4 (3-21) 
t'=2x 10-4 (3.22) 
= 2°6 x 10° years. (3.23) 
This high density appears to predict a very short future for the universe. 
. * (0-01) —2* 10-4 (3.24) 
¢ 2 
=O0°0155 
tg = 0°20 x 10° years. (3-25) 
Here again the time scale is on the small side. The age of the universe for various 
present mean densities is given by the following table 
Po (gm/cc) (t,/T) ty (years) 
° I 13 x 10° 
1:08 x 10°” 0°90 11°7 x 10° 
1:08 x 10° *” 0°67 8-7 x 10° 
1:08 x 10°* 0°30 3°9 x 10° 
1:08 x 107% O°0155 o°2 x 10° 
The range of applicability of these calculations is given by (V 27yp,7r)/c small. 
Taking the value 1-08 x 10°-*° gm/cc for pp, 
27 Po = 3/(477) (3-26) 
implying for applicability that 


(rv 3)/(2zer) <1 (3.27) 


2 
i.e. that r€ — ct (3.28) 
V3 
= 15 x 10° light years. (3.29) 


This distance is considerably greater than the distance of observed nebulae (4) 
and so the approximations involved in this paper are valid. 
4. Suggested further work.—lf 

t=T7/(47ypo)'* (4.1) 
r=CcR/(47yp9)"* (4.2) 
p=poP (4-3) 
v=cV (4.4) 
f=c(47ypo)'*F (4-5) 

equations (1.7) (1.8) (1.10) can be expressed non-dimensionally as 

: @ 2] > ~ 

bare F)=-P+1iF? ry. 

Rap RP) , (4-6) 

or =PV 
Ve eV 


ha i. 
a" a 
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These form a set of three partial differential equations in three variables, and are 
amenable to numerical analysis techniques on the modern large scale computers. 
Basically the interest is in P and V, F being an auxiliary, and it should be 
possible to solve equations (4.6)~(4.8) to provide various model universes. 
Power series solutions of equations (4.6)-(4.8) may be obtained as follows :— 
Let 
P= > P,(T)R™ (4-9) 


n=0 


F= 3 F(T)R™# (4.10) 


ne 
V= > V,(T)R*, (4.11) 
ne QO 


Equations (4.6)-(4.8) become respectively : 


* 


2) 2 @ nel 
Y (2n+3)F,R™+ > P,R™=3 5 DY FF, AR — (4-12) 


n=O n=O nele=l 


> Fy Rett = a a. 


n=0 @=20 


n—V Rent (4-13) 


s F Ret 5 V,,’ Reeth + > (n+1)R***1 s VV, = (4.14) 
n= WU n= n= 0 @=0 
The dash denotes differentiation with respect to 7. The leading terms give 
3Fo= — Po (4-15) 
F,’ = PV, (4.16) 
Fy=V,'+V? (4-17) 


and the other terms may be rearranged to give 


n— 1 
(2n+3)F,+P,=3 > FFi oa (4.18) 


n—1 
F,’-P,Vo-VoPn= 2 Pr-V, (4-19) 


s=l 


n—i 
F,,—V,,' —2(n+1)VoQV, =(n+1) > V,V,_,. (4.20) 
s@=l1 
In equations (4.18)-(4.20), the coefficients of order m are expressed in terms of 
the coefficients of order (n—1), and thus a solution can be obtained. 

Equations (4.15)-(4.17) are the non-dimensional forms of the equations 
dealt with in section 2. One further point of interest arises from equations 
(4.9)-(4.11). The non-linearity of the recession velocity-distance relationship 
is usually expressed 

v=kr+lr?+.... (4.21) 


From equation (4.11), it would appear that a more suitable relationship is 
v=kr+lr+... (4.22) 


and of course, for the integration of equations (4.18)-(4.20) these higher 
coefficients would have to be determined. 
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ANGULAR MEASUREMENTS IN OBSERVATIONAL COSMOLOGY 
W. Davidson 
(Received 1959 August 20) 


Summary 

Theoretical formulae are derived which, by expansion in series, relate the 
observed angular diameters and angular separations of galaxies, radio sources, 
and clusters in an expanding universe to either the measured red shift or 
registered apparent magnitude of representative sources in the region of 
survey. 

It is shown that in principle the cosmological problem may be solved by 
fitting the observed data to these relations. In particular an evolutionary uni- 
verse may be distinguished from a steady-state universe, and the Hubble 
parameter and acceleration parameter for either universe may be determined. 
Allowances are made for possible evolutionary trends in the sources, and it is 
shown how a progressive variation in intrinsic size may be detected. 

Some exact relations between observables are provided for the steady-state 
model. 





1. Introduction..-Yhe apparent sizes of distant galaxies and of clusters of 
galaxies, as indicated by their suitably defined angular diameters, and also the 
angular separations of these objects, provide in principle very important criteria 
for the resolution of the cosmological problem, and in particular for distinguishing 


between an evolutionary and a steady-state universe. It is therefore very much 
to be hoped that statistical measurements of this nature will contribute reliable 
data towards the solution of the problem. 

In this paper formulae are derived to provide an analytic basis for such data. 
In these relations the angular diameters and the angular separations of galaxies, 
clusters of galaxies, and radio sources, allowing for possible evolutionary trends, 
are connected suitably to the red shift exhibited by representative sources in 
the region of survey and also, alternatively, to the apparent magnitude of these 
sources. 

As in the earlier papers in this series on observable relations in cosmology, 
our formulae are designed to connect observables that are, in principle, directly 
measured, and to take account of possible systematic changes with epoch. The 
previous papers, which we shall refer to as I (1) and II (2), dealt with the red 
shift—apparent magnitude relation for galaxies (1), and with the number counts 
of galaxies and radio sources (II). ‘The present paper caters for an observational 
programme that may be regarded as either auxiliary to the programmes covered 
by these papers, or as an independent approach to the solution of the cosmological 
problem. 

It should be mentioned that a valuable paper, dealing specifically with the 
ratio of the angular separation of neighbouring clusters of galaxies to their angular 
diameter as a function of red shift, has been published recently by Florides and 
McCrea (3). We shall have occasion to refer to this paper in our own work here. 
We attempt in this paper, however, to cover a wider programme of observational 
tests, including those that are independent of the requirement by Florides and 
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McCrea that the objects surveyed should not be gravitationally bound to one 
another. Also, it will be a feature of this paper that the intrinsic size of an object 
is regarded as in general a function of epoch, and we shall in principle be able to 
discover by our observable relations whether in fact a dependence on epoch 
exists. In addition, as well as distinguishing between an evolutionary and a 
steady-state universe, the relations found will allow us to identify certain 
characteristic parameters of the particular model that best fits the data. The 
test of Florides and McCrea does not distinguish between the various evolutionary 
models; this was, of course, one of the advantages claimed for their test when it 
was originally suggested by McCrea (4), namely that it is independent of a 
knowledge of the geometry of space time. However, if the universe is evolutionary 
a knowledge of its particular character will obviously be required, and the present 
paper goes some way in providing a means for its determination. 

For the reasons stated in | and II we shall consider it adequate for our purpose 
to adopt the well-known Robertson—Walker space-time metric of an expanding 
universe, in which space is homogeneous and isotropic and a cosmic time exists, 
in the form 
R*(t)(dr? + 7? dé? + r* sin? 0 dp*) | 

(r+ hr®/4)* (1.1) 
Here k=1, 0, or —1 according as the uniform curvature of space (t= constant) 
is positive, zero, or negative respectively. 

According to the usual interpretation of this ideal universe its material contents 
move on the geodesics r=const., 9=const., 6=const. In order to correspond 
with the observed phenomena, however, we shall suppose that, while the centres 
of clusters of galaxies move on these geodesics and thus take part in the expansion 
of the universe, the individual member galaxies will not necessarily do so, but may 
be gravitationally bound within the cluster. ‘The 7, 6, ¢ coordinates of such 
galaxies would in this case vary with the epoch ¢. 

It will be convenient, as in I and I], to change the metric (1.1) to the form 
R(t) f dp” + p* seat sin? 6 dd? | 

Rel asta? (1-2) 





dS? = c* dt? — 





dS? = c? dt* - 


by means of the substitution 
p= Ry, (1.3) 

where R, = R(t) the value of the function R(t) at the present epoch of observation 
t, (fixed). We see that the constant a? is equal to 4R,?/k. 

We repeat here the series expansions of relevant cosmological parameters 
which were established in I: 

If t is the epoch of emission of the radiation from a distant source received at 
the epoch of observation f,, we may write 


R(t)=R,(1 — «4,7 + 4a,.7? — da, (1.4) 


where r=t,—?t and 2,=R,/Ry, %=R,/Ry, etc., suffix o denoting quantities 
evaluated at the epoch ¢,, and a dot indicating differentiation with respect to ?. 

In this notation we may write for the radial coordinate p of the source at the 
epoch of emission 


p=c [1+ a r+ {aay %_) + eI we PETA (1.5) 
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The red shift ratio 5 (=(A,—A)/A) of the radiation from a distant source, 
arising from the expansion of the universe, is given by 
a a oe 
R(t) 
On inverting the series expansion of 5 in terms of 7 we obtain for r in terms of 5 


on (=) 5+ (2) 82+ O(5%). (1.7) 


ty 1 


é I. (1.6) 


This will enable us to eliminate the parameter 7 in favour of the observable 5 
when necessary. 

2. Angular diameter related to proper diameter.—Consider in the first instance 
a distant radiating object which is ideally spherical according to local observers 
on it. Let its proper diameter according to such observers be d and let the radial 
coordinate of its centre be p, both these quantities applying to.the epoch ¢ of 
emission of the radiation that is registered by the observer at p=o at epoch ft. 
We do not assume that d and p are necessarily constant for the sphere. Let the 
angular diameter of the sphere as registered at p=o be @. Since our ultimate 
interest will be in galaxies and clusters of galaxies, or radio sources, whose 
intrinsic dimensions will be small compared with their distance from the observer, 
we shall regard @ as a small quantity of the first order. 

Referring now to the metric (1.2) we may write to a first order approximation 
adequate for our purpose 


re R(t)pé 
Ro(1 + p?/a?)" 


(2.1) 


By this expression the proper diameter of the sphere at the instant of emission 
is related to its angular diameter registered by the observer. 


Discussion 


(i) Application to galaxies or radio sources.—It should be noted that the 
quantity d is not directly measurable by our observer at p=o, but d will have a 
statistical interpretation as a parameter of the cosmological model that we are 
trying to identify. 

In the case of galaxies, because of their varying intrinsic size and the fact 
that they will present varying aspects to the observer, the proper ‘diameter’ d 
requires definition. In any observational test it is to be envisaged that a particular 
class of galaxies, elliptic say, would be selected. As a simple working procedure 
one might perform the angular measure in the direction of the greatest extension 
of the galaxy, as seen in projection, out to that distance on either side of the centre 
where the brightness fell to a certain fraction of that at the centre. There would 
thus be a mean value 8 for the angular ‘diameter’ of all galaxies of the class 
in a given neighbourhood, i.e. within certain limits of exhibited red shift or 
apparent magnitude. There would then be a corresponding mean value d for 
the proper extension of these galaxies at right angles to the line of sight. Ona 
statistical basis this quantity would give a satisfactory measure of the intrinsic 
size of the average galaxy of the specified class at the epoch of emission. 

For galaxies, therefore, our formula will be written in the form 


R(t)p9 


= Ror + pas)’ (2.2) 
19 





274 W. Davidson Vol. 120 


with the above interpretation of d and 6. Ina restricted neighbourhood, at 
sufficient distance, both ¢ and p will vary negligibly from galaxy to galaxy. For 
instance, if the neighbourhood be a whole cluster of galaxies we may with sufficient 
approximation take ¢ and p to correspond to a galaxy at the centre of the 
cluster. 

In the case of galaxies of the specified class selected from clusters selection effect 
could be minimized by restricting the choice to the mth brightest members of the 
cluster, where n is suitably chosen. 

It is clear that in principle all this procedure could be applied to radio sources 
also. It seems probable, however, that most radio sources are too distant for their 
angular diameters to be measurable. It may also be difficult to detect clusters of 
radio sources for the purpose of applying our procedure against selection effect. 
Perhaps the most feasible test with radio sources will be to measure their angular 
separations as a function of their apparent magnitude (Section 5 (ii), (iv)). 

It is important to observe that if an adequate sample of sources were chosen 
for measurement at all distances so that selection effect is eliminated, d could 
be expected to be independent of distance in the case of a steady-state universe. 
Thus, any systematic dependence of d on distance inferred from the observational 
data would confirm that the universe was evolutionary. 

It is of course recognized that angular measurements applied to individual 
galaxies may present formidable difficulties at the great distances at which a 
systematic variation in intrinsic size could be expected to be detectable, and 
this seems certainly so in the case of radio sources. However, it must be 
emphasized that in the case of individual sources our formula is not primarily 
designed for this purpose, but rather to indicate a definite manner of variation 
of 8 with corresponding red shift 5 (or apparent magnitude), via the parameters 
R(t) and p, depending on the cosmological model adopted (Sections 4, 6). It may 
be possible to acquire reliable data for this purpose, from optical surveys at least, 
free of selection effect. Failing this, we may be able to determine the same 
information by studying the dimensions of the clusters themselves, as described 
below, or by the other observational tests introduced in Section 3. 

(ii) Application to clusters—To form a criterion for the intrinsic size of a 
cluster of galaxies (and possibly of radio sources) one might consider all the 
nth brightest members, 1oth brightest say, of the cluster. ‘These could be 
expected to be grouped randomly about the centre of the cluster. As a simple 
definition of the angular diameter of the cluster one might take the mean 6 of 
the angular separations of the five most widely separated pairs of the group, no 
single galaxy being used more than once for this purpose. There would then 
be a corresponding mean proper distance d separating these five pairs, measured 
at right angles to the line of sight, which to adequate approximation would be 
related to 8 in accordance with (2.2). In that formula t and p would again be 
associated with the centre of the cluster to this approximation. In accordance 
with the hypothesis stated in Section 1 the p coordinate of the centre of the 
cluster remains constant, independent of epoch. 

The quantity d, thus defined, could be taken to be a suitable statistical measure 
of the intrinsic size of the cluster. One would naturally choose n fairly large 
for statistical reasons, but small enough to permit the nth brightest galaxies of 
a statistically sufficient number of clusters in each neighbourhood of clusters to 
be observable. Once again, any systematic variation of d inferred from 
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observational data would point to an evolutionary universe. In addition, a 
definite cosmological model must exhibit a definite dependence of 9 on red shift 
or apparent magnitude. 

3. Angular separation of neighbouring objects related to their proper distance 
apart.—Another important criterion for the cosmological problem besides the 
apparent size of the objects is their angular separation in space. As mentioned in 
Section 1 the angular separation of clusters has been treated by Florides and 
McCrea. In deriving a specific formula these authors took the view that any 
variation in the intrinsic size of the clusters with epoch could be justifiably 
neglected. In any case, they argued, such a variation would not in general prevent 
a distinction between an evolutionary universe and a steady-state one. 

In the present paper, however, as in I and II, explicit allowance is made for 
possible systematic evolutionary trends, since there appears every likelihood that 
such important criteria will ultimately be observable. As an example we mention 
the recent exploratory analysis of the population of clusters by Just (5), which 
might be interpreted as tentative evidence for a statistical dependence of cluster 
size on epoch. 

The size of an object has therefore been regarded as a function of epoch in 
general, and our working definition of intrinsic cluster size in Section 2 (ii) 
should reliably take account of any evolutionary variation. We shall also want 
to consider the separation of galaxies (or radio sources) which are gravitationally 
bound within the clusters as well as that of the clusters themselves, and the 
following analysis, although qualitatively similar to that of Florides and McCrea, 
will be rather different in detail. 

Consider a radiating object A, whose radial coordinate is p at the epoch ¢ of 
emission of the radiation that is registered by the observer O (p=0) at epoch fy. 
Let there be another similar radiating object B which is A’s nearest neighbour 
of the same class in the direction making angle ¢, as measured by A, with OA. 
Let the angular separation of A and B as seen by O be ©. Suppose also that the 
proper distance AB at the time of emission is D. (This symbol represented 
‘luminosity distance’ in I but there need be no confusion here.) We shall 
regard D as small compared with OA and hence © as small of the first order. 

Reference to the metric (1.2) then yields, to a first order approximation 
adequate for our purpose, 


R(t)pO 


Dene Ba. 


(3-1) 

We now wish to find the mean value of ©, supposing that D obeys a certain - 
distribution law and the direction of B from A is arbitrary, all directions being 
equally likely. We first suppose that D is constant but ¢ and the plane OAB 
are arbitrary. Clearly, the mean value of © in this case will be 


4a 
© mean = a | O dw, (3.2) 
47 Jo 


where w is the solid angle 27(1—cos¢). ‘Thus 


R(t + p?/a?) aD 
~— Cs” (3-3) 


If D is now supposed to follow its statistical distribution law, of mean D, we 


Omean = 


20 
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find that the mean angular separation of A from its neighbours at epoch ¢ is 0, 
given by 

aD R(t)p© 

4 Rot pla)’ ie 

It is to be noted that while din Section 2 is a mean projected length measured 

at right angles to the line of sight, D on the other hand is the actual mean distance 
separating neighbouring objects. ‘That is, in this case 7D/4 is the mean projected 
distance at right angles to the line of sight. 


Discussion 


(i) Application to individual galaxies or radio sources.—-The substitution of a 
galaxy (or radio source) for object A, and the neighbours of this galaxy taken 
in turn for object B, requires care to guard against selection effect. It would be 
necessary, once again, to confine the examination to the nth brightest members 
of a cluster, with n suitably chosen and taking the galaxy A to be a central member 
of the group. Clusters of radio sources may not of course be detectable as such, 
but perhaps the problem of selection is not so acute in this case since radio 
sources are far less numerous than galaxies. 

‘The manner of the dependence of © on red shift or apparent magnitude, via 
the parameters R(t) and p, is dealt with in Sections 5, 6, and in principle yields 
another method of identifying a cosmological model for the universe. It may 
also be possible to detect a progressive variation in the parameter D, even although 
the sources may be gravitationally bound within the cluster. ‘This of course 
would be incompatible with a steady-state universe. 

(ii) Application to clusters.—In the case of clusters we require the mean 
value © of the angular distances between the centre of a distant cluster A and the 
centres of its immediate neighbours. As a practical procedure one would define 
each cluster by its mth brightest members, for some n. For the angular distance 
between clusters A and B one would take the mean of the angular distances 
between each member of A and a corresponding member of B. 

The parameters ¢ and p in (3.4) will be associated with the centre of cluster A. 
According to the hypothesis in Section 1 the cluster centres take part in the expan- 
sion of the universe, the centres having constant coordinate p. It follows from 
(1.2) that, in an evolutionary universe in which clusters are conserved, the mean 
distance D between neighbouring clusters at epoch ¢ will bear to their mean 
distance D, at epoch f, (the present) the ratio 


(3-5) 


Thus, for clusters in a given distant neighbourhood of an evolutionary universe, 
(3.4) may be written 


nD, p& 
+ I+ p* a* ‘ 


(3.6) 


where p is the radial coordinate of the centre of the neighbourhood, and D, is 
the mean distance separating clusters at epoch ¢,—in particular those in the 
observer's own neighbourhood. 

In a steady-state universe, on the other hand, the mean distance between 
clusters must remain statistically constant for all epochs, due to the steady formation 
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of new clusters in the space between receding clusters. ‘Thus, in a steady-state 
universe, 
D = D, = constant. (3-7) 
These results will form an important criterion for distinguishing between an 
evolutionary and a steady-state universe (Section 5). 


4. Observational tests involving angular diameters 

(i) Angular diameter related to red shift.—The formula relating the proper 
linear diameter (at the epoch of emission) to the observed angular diameter of 
galaxies, radio sources, and clusters, is (2.2), the interpretation of the symbols 
being given in Section 2 (i) and (ii) as appropriate. 

We may now expand R(t) and p in powers of the parameter 7 according to 
(1.4) and (1.5), and in turn in terms of the red shift 6 by means of (1.7). We 
thus obtain 


logy, (50) loge (2) +0°217 (3 ~ 2s) t.0(6?). (4.1) 
1 


The terms of order 5* are not specified because of their complication and the 
fact that they introduce two mofe unknowns whose evaluation by means of the 
observational data would be impossible. However, we are at liberty to formally 
allow for such a term in the analysis of the data by a least squares fitting, as des- 
cribed below. 

The quantity «, is the important acceleration parameter of the expanding 
universe, «, being the Hubble parameter (1). It is evident that if d were statis- 
tically constant, independent of epoch or at any rate varying too slowly to be 
detectable, then the graph of logy (50) against 6 would, to a first approximation, 
be a straight line of gradient 0:217(3—,/%,"). In this case the value of the 
dimensionless parameter ,/«,2 would be immediately determined from the 
data fitted by least squares to a relation of the form 


log, (50) = A + BS + C8*. (4.2) 


If x,/x,” turned out to be sufficiently different from unity (the value for a steady- 
state universe) then an evolutionary universe would be established. 

It might be, however, that d varied sufficiently fast with epoch so that when 
expanded in series it contributed a term of significant magnitude to our first 
order term in 6. Thus, putting 

d=d,- d,t + O(7*), 


where 


ek 


| dt 


d - {40 ; 


and using (1.7), we find 
login (2) = logio (2°) ~ 0°434 (3) 5 + O(8*). (4.3) 
¢ ¢ 14d, 


Here d, is the mean proper diameter for the sources at the present epoch fy. The 
relation (4.1) now becomes 


log, (58) =logyo (24) + 0°217 (3 = “$ ie 2d, 5 + O(8?). (4-4) 


| Hy 
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By (4.2) and (4.4) the data would yield values for 


A=log,, (22) (4.5) 


B=0°217 (3 a oe 2) ; (4.6) 
. 1d, 

By (4.5) a knowledge of either the Hubble parameter x, or of d, would allow 
a value of the other to be deduced, a result which might be of some interest. 

By (4.6) we see that, although we might definitely eliminate the possibility of 
a steady-state universe by studying the gradient of the curve of log (80) against 8, 
we could not establish it with certainty. The data might represent an evolutionary 
universe whose parameter «,/x,? we could not be sure of because of the term in 
d,/d,. Fortunately, however, an independent means of deriving 2,/x," is provided 
by studying the angular separations of clusters, as shown in Section 5. The 
data on angular diameters will then lead via (4.6) to a value for d, d,, assuming 
that we know «,; this will have its interpretation in terms of evolutionary trends 
in the galaxies. 

(ii) Angular diameter related to apparent magnitude.—Since in the case of 
radio sources it is not at present possible to obtain their red shift directly, we 
give an alternative procedure involving apparent magnitude which may also 
be used for galaxies. Some preliminary measurements relating the angular 
diameters of galaxies, measured photoelectrically, to their apparent magnitudes 
have been reported by Baum (6). 

The relation between the red shift 5 and the apparent magnitude m of a 
distant source was derived in I and takes the form 


m= My— 5 logye %4— § + 5 logje (c5) + 1-086 (1 +2 +He+ne) 34 0(82). (4.7) 
1 
Here m is the apparent magnitude as directly registered (apart from aperture or 
obscuration corrections). , is the registered absolute magnitude of a standard 
local source of the same type. ‘The term in «* represents the linearised red shift 
‘correction’ which may be derived by a study of the standard spectrum. ‘The 
term in A* allows for evolutionary change in the spectrum of the distant source 
compared with the standard source. 
On eliminating 6 between (4.4) and (4.7) we find the relation between 6 
and m to be 


log, 9 = — 1 +logy)d, —0°2(m— M,) 


ST 4 eae 2d, 0-2m—M,) + Cf po0-4om—M)} 8) 
+ —(4+n*%+Ae— 10 + O}10 , & 


Therefore, if the observational data for 8 and m were fitted by least squares to a 
relation of the form 


log,)9=E- login ¥ + Gy + Hy?, (4.9) 
with y= 10°%™"-™), 17, assumed known, we should obtain values for the constants 
E=—1+logyody, (4.10) 


os 2178 (44 04n0— 20) inst 
€ ad, 
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The relation (4.10) would give a value for the mean diameter d, of the sources, 
as defined in Section 2, at the present epoch. As regards the value for (4.11), 
we observe that in a steady-state universe, supposing that selection effect had been 
eliminated, we should have A*=o and d,=o. Therefore, given the values of x, 
and «*, we know what G should be for a steady-state universe. If G turned out to 
have a significantly different value then the steady-state model would be ruled out 
in favour of an evolutionary model. On the other hand a steady-state universe 
could not be definitely established on this test alone. 

The disadvantages of using this method for galaxies are of course the fact 
that the parameters .V/, and «*, as well as «,, have to be known, and that it does not 
distinguish between different evolutionary models. On the other hand apparent 
magnitudes are more easily determined than red shift. 


5. Observational tests involving angular separation 

(i) Angular separation of galaxies related to red shift.—The observed mean 
angular separation between a galaxy A and its neighbours is related to its properly 
measured mean distance from them at the epoch of emission by (3.4). This 
relation is formally the same as for angular diameters except for the extra constant 
factor on the left. Analogously to (4.4), therefore, we derive the following 


equation 
logy (50) = logig (~e) +0°217 (3 _ | _ 2D, ) 5+ O(8?). (5.1) 
4¢ 4 


ay Ky ly 


Here 4 is the red shift of the galaxy 4, D, is the mean proper distance separating it 
from its neighbours at the present epoch ft), while Dd, is the present rate of increase 


of this distance with epoch (D, may of course be negative). 
The © —6 observational data would thus provide values for the constants 


(5.2) 


a). (5.3) 


x,Dy 

Eqn. (5.2) allows us to deduce D, if we know «, and vice versa. In a steady-state 
universe D, will be zero statistically, in the sense that galaxies of a given class must 
always be found at the same average distance apart if the average for a region 
(epoch) be taken over a sufficiently large number of clusters belonging to that 
region (epoch). Consequently, since «,/x,? is unity for this case, a coefficient of 5 
found from the data to be significantly different from 0-434 would rule out the 
steady-state model, and establish that the universe was evolutionary. 

In an evolutionary universe ), may or may not be zero, even if the galaxies are 
gravitationally bound in clusters. The assumption that D, is zero would permit 
the value of the important parameter ,/x,? to be deduced. Once again, if we do 
not make this assumption, then a steady-state universe could not be definitely 
established and the value of «,/x,? must be derived from the corresponding data 
on the separation of clusters, as described in (iii). Substitution of this value for 
a,/%,7 into (5.3) would then provide a value for D,/D,, given a, This would 
represent valuable information as to evolutionary trends within the clusters. 

(ii) Angular separation of galaxies or radio sources related to apparent 
magnitude. —To cater specially for radio sources and as an alternative treatment for 
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galaxies we may derive a relation, analogous to (4.8), connecting © to apparent 
magnitude m, in the form 
7D, 


logy) 9 = — 1 + logy, (=) ~o-2(m—M,) 
4 


, 207%, (4+ he 2) 109 20m —M) 4 Of yo9-Mm—M (5.4) 
c a,b; 

Discriminative comments apply to this relation analogous to those made in the 
case of (4.8), where we now read ‘distance of separation’ for ‘ diameter ’. 

(iii) Angular separation of clusters related to red shift.—In the case of clusters 
we have shown in Section 3 (ii) that, because of our hypothesis that the cluster 
centres have constant reference coordinates p, 9, ¢ and thus take part in the 
expansion of the universe, the formula (3.4) is consequently qualified by (3.5), 
leading to (3.6), in the case of an evolutionary universe, whereas for a steady-state 
universe (3.4) is qualified only by (3.7). 

Employing the expansions (1.5) and (1.7) we therefore find that in an evolu- 
tionary universe the mean angular separation ©, defined as in Section 3 (ii), of a 
cluster A from its neighbours is related to the red shift exhibited by a central 
galaxy of A by the equation 

logy» (80) = logy, (a) +0°217 (: - 4) 5 + O(8?), (5.5) 
, 2 
where D, would be the mean distance separating A from its neighbours at the 
present epoch fp. 

On the other hand, in a steady-state universe (3.4) combined with (3.7) 
yield us (5.1) again, except that we have now to put D, =o and 2,/4,2=1. Thus, 
for the steady-state model 


; a,D js 
logy) (80) = logy» (=~) + 0°4348 + O(6?), (5.6) 
4 

where D, is the mean distance separating a cluster from its neighbours at the 

present epoch and at any other epoch. 
On fitting the © —4 data for clusters at successive distances from the observer 

to a relation of the form 

logy) (60) = A” + B’6+ C'S (5.7) 

we would find values for the constants 


A’ =logyo (2) (5.8) 
4¢ 


and 


f % 
» | Or217{t— —=) in an evolutionary universe 
B < a,* 7 

| 


| 0°434 in a steady-state universe. (5.9) 

From (5.8) given the Hubble parameter x, we may calculate D, for either 
kind of universe. 

The eqn. (5.9) is of the utmost value since it directly provides the acceleration 
parameter «,/a,”, whose sign determines whether the expansion of the universe 
is speeding up or slowing down at the present epoch. Unless B” turned out to be 
close to 0*434 the universe would be established as evolutionary. It is to be noted, 
however, that even if B” was near to 0-434 it might signify an evolutionary universe 
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for which «,/«,2 was equal to —~1. ‘To remove ambiguity in this case one could 
appeal to the other tests described in this paper, or of course to the tests described 
in I and Il. 

(iv) Angular separation of clusters related to apparent magnitude.—Assuming 
that radio sources occur mainly within the clusters of galaxies, in accordance with 
the collision hypothesis for their principal origin, it may be useful to provide a 
formula for the angular separation of clusters in terms of the apparent magnitude 
of representative radio sources within them. As mentioned in Section 2 (i), 
it may be difficult to identify clusters of radio sources, but if radio sources occur 
only in clusters this would not matter. The formula will in any case be suitable 
as an alternative treatment for clusters of galaxies observed optically. 

In the case of an evolutionary universe it is clear that the required formula will 
simply be (5.4), where D, is now the mean distance between cluster centres at 
the present epoch and, in virtue of (3.5), Dy=%,D). Hence 


log, = — 1+ logy, (=~) —0-2(m— My) 
4 
— ria (2+ «*+A*)10° 2" —Mo) + Of 10° MM}, (5.10) 


For a steady-state universe we have to put D,=0. In addition A* =o, since 
statistically the sources will be of strength independent of epoch. In this case, 
therefore, 


logy = ~1+logiy ("—) ~o'2(m—M,) 
4 


2°17 , 
+ oe (4+ « %) 10° %Xm—Me) 4+ Of 109 Mima), (5.11) 
Cc 


The fitting of the data to a relation of the form 
logy = E’ —logyy + G’y + H'y*, (5.12) 


where y = 10° %™"~™»), would present us with values for the constants 


E'= 1 + logy (=), (5.13) 


f 
2°17% : ‘ ; 

|; £1 (2+«*+A*) inanevolutionary universe 

” c 

G'= (5-14) 


217% ; : 
—__! (4+«®) in a steady-state universe. 
- 


Eqn. (5.13) gives us the mean distance between cluster centres at the present 
epoch. Egn. (5.14) indicates that there is a definite value for G’ in a steady-state 
universe which is ascertainable since «* and «, will be known. Coincidence 
apart, therefore, the two kinds of universe should be readily distinguishable. 
This would be especially so if A* is negative for an evolutionary universe (average 
power of the sources diminishing as the epoch advances). 

(v) The test of Florides and McCrea.—As a measure of the congestion of the 
universe in a given region, Florides and McCrea (3) suggest dividing the mean 
angular separation of the clusters in that region by the mean angular diameter 
of these clusters. Except in a special circumstance this ratio will be a function of 
distance in an evolutionary universe, but constant in a steady-state universe, so 
that this constitutes a valuable means of distinguishing between these two kinds 
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of universe. According to the analysis developed in this paper their test may be 
expressed as follows. 
For clusters in a given neighbourhood of an evolutionary universe we get, 
dividing (3.4) by (2.2), 
oe) " aD 


6 4d 
(Our practical procedure for finding 9 and 0, as defined in Sections 2 and 3, is 
rather different from that of Florides and McCrea but this need not concern us 
here.) Also, the mean distance D separating a cluster from its neighbours at the 
epoch of emission is related in accordance with (3.5) to the mean distance separating 
them at the present epoch t,. Hence 
6) ~ 7D, R(t) 
6 4d : Ry 
It follows that in an evolutionary universe the ratio of the mean angular separation 
of a cluster A from its neighbours to the angular diameter of A may be written, 
observing (1.6), 


(5.16) 


(5-17) 


where 4 is the red shift exhibited by the central region of A. 
On the other hand, for a steady-state universe dividing (3.4) by (2.2) and 
putting D=D,, d=d,, leads to 


(5.18) 


It follows that unless a possible variation of din (5.17) exactly masks the trend 
in the factor 1/(1 +4) the ratio of the observables must vary systematically for an 
evolutionary universe while remaining constant for a steady-state universe. 

Should the ratio turn out to be in fact constant any ambiguity could only be 
settled by appealing to another test. In any case, recourse to other tests would be 
necessary even if the ratio were not constant, if we wanted to distinguish between 
different evolutionary models by determining kinematic parameters such as 


Oty /a,?. 


6. Some exact observable relations for a steady-state universe 


(i) Introduction.—For the steady-state model the metric factor R(t) is equal 
to e'T, where T is a constant which is evidently the reciprocal of the Hubble 
parameter a, for the model. In addition the curvature of space is zero so that 
1/a2=o0. It is therefore possible to give exact expressions for some of our rela- 
tions between observables in this case, always bearing in mind that our formulae 
are subject to the restriction that the mean angular diameter 6, or the mean angular 
separation 0, of the objects is a small quantity, so that an approximation in this 
respect has already been made. 

It will, unfortunately, not be possible to express angular diameter and angular 
separation in terms of apparent magnitude by an exact relation. This is because 
the exact red shift ‘correction’ to registered apparent magnitudes, discussed in I, 
contains 6 explicitly, making necessary a measurement of red shift as well as 
apparent magnitude. If we wished to confine distance measurements to apparent 
magnitude only, we must therefore fit the data to the approximate relations already 
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dealt with. In these, of course, 5 has been eliminated in terms of apparent 
magnitude by a first order approximation. 

Referring to (1-6) we see that in the steady-state model the red shift 5 is related 
to the epoch ¢ of emission of the radiation by the equation 


§=e4-OT_ 1, (6.1) 


where as usual ¢, is the epoch of its reception. 
Also, the null geodesic followed by the radiation between its source, of radial 
coordinate p at the moment of emission, and the observer at p = 0 has the equation, 


in the general case, 
re = dp [" dt 
an oe eT) er 6.2 
| repre : ° J, R(t) i 
In the case of the steady-state model this leads to 
p=cT (e&-%T_ 7), 
or p=cTS, (6.3) 
by (6.1). 
(ii) Angular diameter related to red shift.—By (2.2) and (6.3) we obtain for a 


steady-state universe the exact relation between the mean angular diameter of an 
object and its red shift in the form 


_ d(1+8) 
6= rr. (6.4) 


where we have put d= d, as dis statistically independent of epoch in this model. 


We see that in a steady-state universe the mean angular diameters of galaxies, 
radio sources, and clusters of these objects will decrease monotonically as their 
exhibited red shift increases, towards a definite lower limit d,/cT at the observa- 
tional horizon of the universe (6= 0). This limit is a calculable quantity if we 
ascertain d, by local observations and know the value of the Hubble parameter 1/T. 

The result (6.4) has been derived and discussed by the author previously (7), 
when an estimate was given for the limiting angular diameter of a standard galaxy 
in a steady-state universe. Since then the Hubble constant has been revised, 
the present estimated value for its reciprocal being 13 x 10°yrs. Taking this 
value for T we find that the limiting value of 9, corresponding to d, = 30,000 It. yrs. 
say, would be approximately 0-5 seconds of arc. For clusters, taking for d, the 
rough figure of 5 x 10° lt. yrs., the lower limit would be nearly 80 seconds of arc. 

Formula (6.4) has been discussed more recently by Hoyle (8). He calculates 
that the limiting value of 8 (our notation) for a radio source, of a standard intrinsic 
diameter equal to that of the source in Cygnus, would be approximately 4 seconds 
of arc. For comparison Hoyle has given an interesting discussion of the analogous 
situation in the Einstein-de Sitter universe, a well known evolutionary model. 
He finds that in this model, neglecting any progressive variation of the intrinsic 
diameter, 6 would decrease to a minimum at 6=5/4, and thereafter steadily 
increase. For the standard radio source described above he finds that the 
minimum value of § would be approximately 15 seconds of arc. 

It may be shown, however, that this behaviour is not general for the evolu- 
tionary models, even in the unlikely event that d remains constant. For instance, 
in a Milne universe 9 would decrease monotonically to the value 2dx,/c at the 
observational horizon, «, being the Hubble constant. This limit is incidentally 
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twice that of the steady-state model. The behaviour of 9 in the evolutionary 
models depends on the metric factor R(t), the curvature of space, and the vari- 
ability of d, so that no useful general criterion is available for these models in an 
exact form. 

The most positive conclusion that we can make, therefore, is that since the 
exact form of dependence of 9 on 8 in the steady-state model is likely to be unique, 
the strict behaviour (within statistical error) of the observables in accordance with 
(6-4) would represent powerful evidence in support of the steady-state model. 
On the other hand, any disagreement of the data with this relation would rule out 
such a universe. 

(iii) Angular separation related to red shift.—Formula (3.4) relates the mean 
proper distance separating neighbouring galaxies, radio sources, or clusters, at 
the epoch of emission, to their observed mean angular separation. For the 
steady-state model, putting D= D,, we therefore find 

7Do =pelt-WT@, 
4 
e (1+8) 
Whence Q= a (6.5) 

Thus, in the same way as for angular diameters, the mean angular separations 
of these objects decrease monotonically as the distance increases, to the lower 
limit 7D,/4cT. Once again this is a calculable quantity. For galaxies of a 
certain order of brightness, separated within the clusters by a mean distance of 
10° lt. yrs. say, the limiting angular separation would be approximately 12 seconds 
of arc. For whole clusters separated by a mean distance of 2 x 10’ It. yrs. the 
lower limit of © would be over 4 minutes of arc. 

In the case of the evolutionary models no general statement of an exact nature 
can be made about the mean angular separation of galaxies within the clusters, 
not least because of the possible dependence of their mean distance D on epoch. 
For the clusters themselves, however, it is possible to be more specific. Since we 
have assumed that the cluster centres follow the expansion of the universe, the 
formula (3.6) applies. The author has found that the possible types of behaviour 
of the mean angular separation © of cluster centres may be categorised as follows. 

(a) Curvature of space positive. In this case, if space is spherical so that 
antipodal points are not identified, then it is possible for © to steadily decrease to a 
minimum value of 7D,/4R, (with our usual meanings for the symbols) occurring 
at p=2R,, and thereafter to steadily increase. However, the limiting value of 
© depends on the location of the observational horizon in the model, which in 
turn depends on the function R(t) in (1.2). If space is elliptic then antipodal 
points are identified, © would be monotonic decreasing and the limiting value of 0 
would be > D,/4R, depending on the location of the horizon. 

(6) Curvature of space equal to zero. For these models © must steadily 
decrease to a limiting value which is > according to the location of the horizon, 
whose coordinate p will be < a. 

(c) Curvature of space negative. In this case © must steadily decrease to a 
lower limit >, occurring at the observational horizon of coordinate p < 2Ry. 

It may easily be verified that in the Einstein-de Sitter universe, which belongs 
to category (b), the limit for © is the non zero value 7D x,/8c. This is half the 
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steady-state value. The Milne universe is of category (c) and in this case the 
lower limit for © is zero. 

We may assert, therefore, that the formula (6.5) represents another powerful 
criterion for the steady-state model, especially in its application to the angular 
separation of clusters. In particular, if beyond a certain distance the angular 
separations of neighbouring sources, or clusters of sources, were observed to be 
consistently less than the respective limits for the steady-state model then that 
universe would be eliminated. 

(iv) The ratio of angular separation to angular diameter.—The expression for 
this ratio, namely 


7 
=—, (6.6) 
4d, 
which was derived and discussed in Section 5(v), is another exact result for the 
steady-state model in the sense described in (i). 
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